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The structure of the 1/Nc expansion for the baryon distribution amplitude in QCD is tested using 
quark models. Earlier conjectures about this structure based on the evolution equation and on the 
soft-pion theorem are confirmed by the model analysis. The problem of the calculation of the baryon 
wave function at large Nc is reduced to the analysis of equations of motion for an effective classical 
dynamical system. 
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I. INTRODUCTION 
A. Large Nc'. QCD and models 



The limit of the large number of colors Nc and the 1/Nc expansion Sixe important nonperturbative methods 

in QCD. However, in spite of the long studies of the large- A^c limit, QCD is not solved even in the leading order of 
the 1/Nc expansion. In the applications of the l/Nc expansion to QCD, one usually deals not with the order-by-order 
calculation of the 1/Nc corrections but with the analysis of the structure of the 1/Nc expansion. 

The failure of the dynamical approach to the 1/Nc expansion in QCD led to the redistribution of the effort towards 
various models imitating large-A^c QCD. This approach proved to be rather fruitful. Instead of being trapped into 
the study of model artifacts (which often happens in the phenomenology of strong interactions), the analysis of 
various large- iV^ models has revealed a number of their common properties including the so-called large- iVc spin-flavor 
*n symmetry l3-[il| which is believed to be a symmetry of large- iVc QCD as well. 

A classical example of results, historically first derived in the Skyrme model but actually based on the spin-flavor 
' symmetry of large- iVc QCD, is the ratio of the pion-nucleon (gyrNN) and pion-nucleon-A (g-^NA) constants jl2| 

Qh 

^- El^ = ^[i + oiN-')] (exp.1.48). (1.1) 

' The same ratio 3/2 can be derived in the quark model 13] in the limit Nc — > oo. 

This simple example shows that the study of large- A^c models can be helpful for establishing the properties of 
large- A^c QCD itself. 



B. Baryon wave function in the large-TVc limit 



This paper is devoted to the 1/Nc expansion for the baryon wave function. This problem was raised in Refs. p^ll5l| . 
The construction of the 1/Nc expansion for the baryon wave function is rather nontrivial, since one has to deal with 
the exponential large- A'c behavior. In Refs. 0,^3 was shown that this exponential behavior can be consistently 
described using a specially designed generating functional for the baryon wave function. Strictly speaking, the results 
of Refs. rpjllTll are based on conjectures, whose consistency was checked using various indirect methods like evolution 
equation jig and soft-pion theorem 17]. In this paper we provide additional arguments supporting these conjectures, 
using models of large- iVc QCD. 

Before describing the main features of the 1/Nc expansion for the baryon wave function we would like to agree 
about the definition of the baryon wave function. The concept of the wave function is natural for nonrelativistic 
quantum mechanics but its incorporation in QCD requires some extra specifications. Usually one keeps in mind the 
transition matrix element between baryon \B) and vacuum ]0) 

{0\X{Mzc)\B) , (1.2) 

c=l 

where ijjc{z) is the quark field with the color index c. Strictly speaking one has to settle several questions: 
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1) ordering of quark operators ipdzc), 

2) insertion of Wilson lines providing the gauge invariance, 

3) choice of the space-time or momentum kinematics (equal-time wave function, light-cone wave function etc.). 
From the point of view of QCD applications to hard exclusive processes the most interesting object is the baryon 

distribution amplitude |1^ ,2^ which corresponds to the choice 

Zc = XcU, = (1-3) 

in Eq. I|1.2(l . Here n is a fixed light-cone vector. The standard definition of the baryon distribution amplitude also 
involves the Fourier transformation from parameters Ac to variables Xc corresponding to the momentum fractions of 
quarks in the infinite momentum frame: 

. AT, 

/ dXe'^'''"'{0\Y[^cd^kn)\B) ^ ec,c,...c^^'fB{xi,X2,...,XNj. (1.4) 
k=i 

We have omitted various normalization factors which are not important for the discussion of the general properties 
of the 1/Nc expansion. 

Although the light-cone elements of the definition of the baryon distribution amplitude are important from the 
dynamical point of view, they are rather insignificant for the structure of the 1/Nc expansion. What is really important 
for the construction of the 1/Nc expansion is that the baryon wave function depends on Nc quark variables. At this 
moment one meets a serious problem: how can one construct the expansion for the function (|1.4|l which depends 
on Nc variables Xfc? This problem is much more general than the context of the baryon distribution amplitude. Below 
we will consider the general case without specifying the precise meaning of the wave function ^> b{xi^ X2, . . . ,xn^) 
and of its arguments xi, X2, ■ ■ ■ , xn^^- For example, the variables Xk can be composite objects carrying information 
about the spin and fiavor of quarks in addition to their coordinates or momenta. Since the antisymmetric color 
tensor Scj^c2...cn is factored out in Eq. (|1.4|) . the function si^i, X2, ■ ■ ■ , xn^) is symmetric in all its arguments 
Xi,X2, ■ ■ ■ ,a;Af ■ 

In Refs. [la. Il7j the problem of the construction the 1/Nc expansion for function depending on Nc variables was 
solved by introducing a generating functional for the baryon wave function. Omitting irrelevant technical details (spin 
and flavor structure, light-cone kinematics, 1/Nc rescaling of the quark momentum, etc.) and concentrating on the 
general structure which is important for the large- A'^c limit, we can define this generating functional as 

*s(g) = j '^b{xi,X2, ■ ■ . ,XN,)g{xi)g{x2) ■ . .g{xNjdxidx2 ■ ■ - dxN, , (1-5) 

where g{x) is an arbitrary "source" function. Since function b{xi,X2, . ■ ■ ,xm^) is symmetric in permutations of 
xi,X2, ■ ■ ■ ,xn^, the transition from the wave function b(xi,X2-, ■ ■ ■ ,xn^) to the functional $5(5) does not lead to 
any loss of information at finite Nc- 

As was argued in Refs. ^^e functional $5(5) has an exponential behavior at large Nc- 



<^B{g) ^^°° e^v[NcW{g)] . (1.6) 

where W{g) is some [Nc independent) functional of g. An important property of the functional W{g) is its universality: 
all low-lying baryons [with 0{N'^^) or 0{N^) excitation energy] are described by the same functional W{g). 
The dependence on the type of the baryon B appears in the preexponential factor AB{g)'- 

^B{g) K^AB{g) exp [NcW{g)] [l + 0{Nc')\ . (1.7) 

Although the functional AB{g) depends on the baryon _B, this dependence has simple factorization properties studied 
in Ref. 

Strictly speaking, these properties of functional W{g) and AB{g) were not properly derived in Refs. 0,^3- 
Instead the consistency of these properties with a number of basic features of QCD was checked: 

- QCD evolution equation [T^ . 

- asymptotic behavior in terms of anomalous dimensions of leading twist operators [T^ , 

- spin-flavor symmetry of large- A^^c QCD j 

- soft-pion theorem for the baryon distribution amplitude . 

Although all these consistency checks successfully passed, they cannot be considered as a final proof. In this paper 
we turn to the analysis of large- A"c models in order to find additional arguments supporting the large- A^c behavior 
(|1.7|) and illustrating the basic properties of functionals W{g) and As(g)- 
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C. Consequences of the spin-flavor symmetry 

In the absence of the dynamical solution for large- A'c QCD, one cannot expect to find much more than relations like 
the A^-A ratio in the pure large-iV^ approach. In Ref. [I3 a relation was derived for the generating functionals 
and describing the baryon wave functions of nucleon and Delta resonance where and J3 are 

projections of isospin and spin, respectively. 

In order to write this relation let us first notice that $^^^(3) is a 2 x 2 matrix with indices T3, J3. Let us consider 
the 2x2 matrix 



R = 



Vdet$W(5) 



(1.8) 



Obviously 



det i? = 1 



(1.9) 



so that matrix R belong to the group SL{2, C). The irreducible representations of SL{2, C) are described by Wigner 
functions D^^^'^'^\R) depending on two "spins" ji, j2- Let us consider the Wigner function corresponding to ji — 3/2, 
22 = 0. For brevity we will denote it simply Dr^l^j^. Now we are ready to write the 7V-A relation derived in Ref. (17| : 



Vdet<i>(^)(5) 



2det$W(5) [l + 0{N-^)\ 



(1.10) 



Relation (fTTIHl was derived in Ref. |17| using the realization of the large- A^c spin-flavor symmetry in the space of the 
preexponential functionals Asig) appearing in Eq. (|1.7|) . 

Due to the universality of the functional W{g), relation H1.10(l is equivalent to the analogous relation for the 

preexponential functionals A^j^((7) and ^^^^(5) appearing in Eq. 11.7|l : 



A(^)(£ 



VdetylW(.g) 



^2detAW(.g) [l + 0(7V-i)] 



(1.11) 



D. What do we want of models? 



It should be stressed that the construction of viable realistic models of the baryon distribution amplitude is not the 
aim of this paper. We are interested in another problem. We want to trace how the exponential large- TVc behavior 
(|1.7|) appears within a certain class of models. The main problem of large- iVc QCD is that its dynamics is not solved. 
Therefore the analysis of Refs. [Til IT^ was based on the assumption about the exponential large- iVc behavior \1.7\ 
with the universal functional W{g). 

In QCD we cannot compute the functional W{g). On the contrary, in the models studied in this paper both the 
exponential behavior (|1.7|l and the calculation of W{g) are under good theoretical control. From the point of view of 
the aim of this paper (model check of the assumptions made in large- iVc QCD), the result of this work is trivial: all 
assumptions involved in the theoretical construction |l(lll7| are confirmed in the models studied in this paper. This 
puts the results of Refs. on more sofid ground. 

Several comments must be made about the status of the spin-flavor symmetry in our model analysis. Most of this 
work is based on the Hamiltonian (Schrodinger) analysis of models. As is well known, the mass difference Ma — Mn 
between the A resonance and nucleon is 0(A^~^), whereas the masses and Ma grow as 0{Nc). Therefore the 
mass splitting A/a — Mn is an effect of the next-to-next-to- leading order of the 1/Nc expansion. In this paper we 
do not go so deeply into the l/Nc expansion. Actually the main part is devoted to the calculation of the functional 
W{g) corresponding to the leading order. Therefore the subtleties of the realization of the spin-flavor symmetry are 
not important for the most of our analysis. Moreover, in our analysis we make no special assumptions about the 
symmetries of the Hamiltonian [only the SU{Nc) color symmetry is important for us]. 

On the other hand, the readers interested in model implementations of the spin-flavor symmetry can be addressed 
to Ref. [T^ where functionals W{g) and Asig) are computed in the naive quark model and the results explicitly 
agree with all constraints imposed by the spin-flavor symmetry. 
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E. Mean field approximation and large-A^c limit 



In contrast to "unsolvable" large- A'c QCD, many large- A^c models allow for a straightforward construction of the 
1/Nc expansion. Although the precise form of the l/iV^ expansion can be model dependent, many features are 
common since they are based on the mean field approximation which is justified in these models (but not in QCD) 
in the large- A'^c limit. The manifestations of the underlying mean field approximation can be different: saddle point 
approximation in the path integral approach, the Hartree equation in the Hamiltonian approach to the models with 
explicit quark degrees of freedom, classical equations of motion in models based on meson fields, etc. 

In this paper we work with models explicitly containing quark degrees of freedom. In the leading order of the 1 /N,, 
expansion the solution of these models is described by the Hartree equation. 



F. Density matrix or wave function? 

This paper deals with the 1/Nc expansion for the functional $B(ff) H1.5|l describing the baryon wave function. 
One can wonder how this approach is related to the mean field method. Note that the mean-field approach to 
fermion systems is based on equations for the density matrix whereas the functional $5(5) is constructed in terms of 
the wave function b {xi , X2 , . . . , a;jVc ) • Strictly speaking the transition from the multiparticle wave function to the 
single-particle density matrix 

p{x\,xi) = j dx2.. .dxN,^B{x[,X2, ■ ■ . ,XnJ'^*b{xi,X2, ■ ■ ■,xnJ (1-12) 

leads to a loss of information. Therefore the knowledge of the mean field solution in terms of the density matrix 
p{x'^, xi) is not sufficient for the calculation of the functionals W{g) and Asig) describing the large- A^c behavior H1.7|) 
of the baryon wave function. Nevertheless one can derive a closed equation for the functional W{g) [See Sec. II VI and 
Eq. H4.6|) ]. An essential part of this paper is devoted to the analysis of this equation and to the construction of its 
solutions. 

In the leading order of the l/Nc expansion the density matrix p(x',x) is described by the Hartree equation which 
leads to the representation for p{x' , x) in terms of single particle wave functions ijjn{x) of occupied states: 

p^^'^^)^^^ ^ V«(^')C(^)- (1.13) 

n:occupicd 

In the simplest models with only one "valence" level V'o occupied by Nc quarks, the density matrix factorizes into the 
product 

pix'^x^^'^^M^'Woix). (1.14) 

Here we assume that the trivial color term Sec' is factored out. 

In the large- A^c limit the factorization of the density matrix H1.14|l becomes asymptotically exact. Naively one could 
think that the wave function \E'(a;i,a;2, . . . ,xn^) has a similar factorization at large Nc'. 

^{xi,X2,.-.,xnJ i/'o(2^fc) (wrong!). (1-15) 

fc=i 

Unfortunately this statement is wrong. In order to see the violation of Eq. H1.15|l let us multiply it by JIa^i 5(^fc) 
and integrate over all Xk- If Eq. I|1.15|) were correct, then we would obtain using Eq. (|1.7|) 



W{g) = In 



ipo{x)g{x)dx 



(wrong!). (1-16) 



In Sec. I VIII we describe a model for which the functional W{g) can be computed analytically [Eq. (|7.74l) ]. This result 
explicitly shows the breakdown of the naive expectation l|1.16|l . 

Although the naive statement (|1.15l) about the factorization of the wave function can be often met in literature, 
this statement can be trusted only in the sense of the factorization of the density matrix H1.14|l in the leading order 
of the 1/A"c expansion. 

To summarize, we cannot use Eq. (|1.16|) for the calculation of the functional W{g). Moreover, the knowledge of 
the solution of the Hartree equation is not sufficient for the calculation of W{g). In this paper we derive a differential 
equation for W{g) and construct its solution. 
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G. Factorization of the preexponential terms 

The universality of the large- A'c behavior (|1.6() holds only with the exponential accuracy. The preexponential 
functional Asig) in Eq. (|1.7|) depends on the baryon (baryon-meson) state B. In Refs. [19. Il7j several important 
factorization properties of these functionals were suggested and checked using the evolution equation, the spin-flavor 
symmetry and the soft-pion theorem. In this paper we verify these factorization properties by explicit calculations in 
large- iVc models. 

One should keep in mind that the factorization properties of functionals Asig) are determined by two ingredients: 

1) the breakdown of the spin and flavor symmetries by the mean field solution and the restoration of these symmetries 
via the standard method of "the quantization of collective coordinates" , 

2) the spectrum of the harmonic fluctuations above the mean field solution. 

If the mean field solution for the baryon had the same symmetries as the vacuum mean field solution, then we 
would have a simple oscillator-like spectrum of baryons 

£b = £{n,} = N,Eo + (^AEo + X^nfcfifcj + 0{N-') , (1.17) 

where the baryon excitations B are labeled by sets of integer numbers {n^} associated with elementary fi^ excitations. 
In this case the factorization of Asig) is described by the formula 

ABig)^A^^^^{g)^A("\g)l[[A,ig)r ■ (1.18) 

k 

In the case when the mean field solution breaks some symmetries (which are nevertheless restored by the collective 
quantization) the situation is more subtle and the precise expression for As (g) depends on the involved symmetries. 
As is well known, the effects of broken symmetries in the mean field approach to large- A^c models correspond to the 
manifestation of the spin-flavor symmetry in large- A^c QCD. The structure of the functional Ab (g) in large- A'c QCD 
and the role of the spin- flavor symmetry were studied in Ref. [T^ . 

In this paper we concentrate on models where the effects of broken symmetries and zero modes are absent so that 
one has the simple factorization formula (|1.18|l . In Sec. IVIII we describe a model where the factorization relation 
(|1.18l) can be checked explicitly. 



H. Schrodinger equation versus path integral 

One can use two methods for the analysis of the models for the baryon wave function: the stationary Schrodinger 
equation or the path integral approach. In principle, these two methods should give equivalent results. It is well 
known that the large- A'c limit has many common features with the semiclassical limit. In the path integral approach 
this leads to a representation for W{g) in terms of classical time- dependent trajectories. In the case of the stationary 
Schrodinger equation one works with the time-independent formalism. 

In the simplest models (i.e. models with the trivial vacuum), the connection between the two representations 
can be easily seen. However, in the general case the equivalence of results obtained in the two approaches is less 
obvious. Therefore in this paper we use both methods. The first part of the paper (Sees. ITTllVIII|l is based on the 
Schrodinger equation. In Sec. IIXI we describe the path integral approach and derive a representation for W{g) in 
terms of trajectories. After that we demonstrate the equivalence of the two representations for W{g) in Sec. |3 



I. Large-A'c limit and classical dynamics 

An essential part of this paper is devoted to the representation of W{g) in terms of classical trajectories described 
by an effective Hamiltonian. In fact, the functional W{g) can be interpreted as a classical action depending on the 
"coordinate" g. This action obeys the classical Hamilton-Jacobi equation. The method of trajectories can be used 
for the construction of the solution of this Hamilton-Jacobi equation. 

The connection between the large- A' systems and classical dynamics was extensively discussed in literature (see 
also review and references therein). In the path integral approach to the large- A" systems, the 1/N expansion 
can be constructed using the saddle point method. In a straightforward approach the corresponding saddle point 
equations have a nonlocal form. Nevertheless in many cases it is possible to find new degrees of freedom which reduce 
the saddle point equations to a local form described by a local Hamiltonian dynamics. 
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In this paper we use similar methods and reduce the problem of the calculation of the functional W{g) to the analysis 
of trajectories described by an effective Hamiltonian. However, the problem studied in this paper is different from the 
"traditional" time-dependent problems in large- A*" systems 21]. Most of the works on the effective classical dynamics 
of large- TVc systems concentrate on the time dependence of expectation values of observables so that the classical 
trajectories correspond to the time-dependent diagonal matrix elements of the corresponding quantum operators. On 
the contrary, the calculation of the baryon wave function is a problem of nondiagonal transition matrix elements. In 
spite of the difference between the two cases, we will see that the corresponding effective classical dynamical systems 
have many common features. 



In our analysis of the large- iVc models we will meet many equations which are well known in the "traditional" many- 
body physics (see e.g. Ref. p^l. The machinery of the old many-body physics includes many interesting equations: 
stationary Hartree-Fock, time-dependent Hartree-Fock equation, random phase approximation (RPA) equations. A 
common feature of these equations is that their derivation is usually based on ad hoc approximations so that in a 
general case the validity of these approximations is a matter of luck. For example, the standard derivation of the 
Hartree-Fock equation is based on the variational principle with a special ansatz for the wave function. Nevertheless in 
some cases one can find a justification for these approximations. In particular, the 1/N expansion (N not necessarily 
being color) puts these equations on solid ground p3. l25l|. In the framework of the 1/iV expansion, the Hartree 
equation corresponds to the leading order of the 1/N expansion whereas the contribution of the Fock term is 1/N 
suppressed. Similarly, in the next-to-leading order of the 1/iV expansion one arrives at the RPA equation |3 EBl 
with slight modifications caused by the large- TV limit 

As was explained above, the functional W{g) cannot be expressed via the solution of the stationary Hartree 
equation. However, in Sec. IVII we will see that W{g) can be expressed via solutions of the so-called time- dependent 
Hartree equation (TDHE). Indeed, as was mentioned above, the functional W{g) is just an action for some effective 
mechanical system. The equations of motion describing these trajectories of this mechanical system are similar to 
the TDHE. The properties and solutions of the TDHE were extensively studied in the context of the problems of the 
many-body physics (typically in its Hartree-Fock version ^2^). Strictly speaking, the analogy between our equations 
and TDHE is not complete since our problem the large- A'c wave functions differs from the traditional problems of 
the many-body physics. Nevertheless many interesting connections can be found. For example, our analysis of the 
effective large- iV^ dynamics makes use of the Hamiltonian structure of the classical equations of motion. In the context 
of the "traditional" approach to TDHE (with Fock term included) , the classical Hamiltonian structure of the TDHE 
was studied in Ref. 1271. 



In the simplest models the baryon consists of Nc quarks without any quark-antiquark pairs. This class of models 
is studied in the first part of the paper (Sees. IIIII - IVIlll . The second part of this paper (Sees. IVIIIKlX|) deals 
more complicated models where the quark-antiquark pairs are described using the well-known "old-fashioned" Dirac 
representation for antiparticles. In these models we have a nontrivial physical vacuum made of MNc quarks put into 
the bare vacuum (with some integer M), whereas the baryon is considered as a state containing (M -I- l)Nc quarks 
above the bare vacuum. The baryon wave function [defined by Eq. (|1.2I) as a transition matrix element between the 
baryon and the physical vacuum] still can be used for the construction of the generating functional ^B{g) H1.5|l . We 
show that in this class of models we also have the exponential behavior (|l.t)|l of ^B{g) and develop methods for the 
calculation of the functional W{g). 



The general structure of the 1/Nc expansion for the functional $5(17) is the same in quantum mechanics and in 
quantum field theory. For simplicity we use the discrete quantum mechanical notation in this paper. For example, 
instead of the continuous variables Xk of Eq. H1.5|l we write discrete indices ik- 



J. Large Nc models and traditional many-body physics 



K. Antiquarks in large-A'c models 



L. From quantum mechanics to quantum field theory 




(1.19) 
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This compact discrete notation allows us to simplify equations. By inertia we will often call $_B(.g) generating 
functional, although the word function could be more appropriate. 

One can ask whether our compact discrete notation is only a matter of language or some serious problems will come 
in the case of field theoretical models. The main problem is that our quantum- mechanical models are based on the 
four-fermion interaction. Their field theory analog of these models is represented by model of the Nambu-Jona-Lasinio 
(NJL) type. Formally our quantum-mechanical equations can be generalized for case of the field theory. However, 
there is a problem with the nonrenormalizability of the four-fermion interaction. From the physical point of view, 
models of the NJL type can be considered only as low-energy effective models which should be simply cut at high 
momenta. But problem is that our treatment is based on the Hamiltonian (or path integral) formalism which assumes 
the locality in time. Therefore applications of our formalism to NJL and similar models would require regularizations 
preserving the locality in time. Unfortunately most of rcgularization local in time conflict with the Lorentz invariance. 

It should be stressed that the methods developed in this paper can be applied only to the large- iVc models of QCD 
but not for QCD itself. Nevertheless some special problems of large- iVc QCD these methods still can be used:: 

1) the problem of heavy-quark baryons 

2) the problem of the diagonalization of the matrix of anomalous dimensions for the leading-twist baryon operators 

M. Structure of the paper 

As was already mentioned, the paper consists of two parts. The first part (Sees. IIIIIIVIip is devoted to models 
with the trivial vacuum (i.e. without antiquarks) whereas in Sees. IVIIII lx| more powerful methods are developed for 
models imitating antiquarks in terms of the old Dirac picture. The general class of models is discussed in Sec. ^ In 
Sec. mil we describe the simplest models of the baryon wave function where the baryon appears as a bound state of 
Nc quarks (without quark-antiquark pairs). In Sec. IIVI we derive a nonlinear differential equation for the functional 
W{g) and show that this equation agrees with the standard Hartree equation for large- systems. In Sec. |3we turn 
to the analysis of the preexponential functional A{g) (|1.7|) . In Sec. I VII we show how the problem of the calculation of 
W{g) can be reduced to the analysis of trajectories in some effective classical dynamics. In Sec. I VIII we study a simple 
model where the functional W{g) can be computed analytically. In Sec. IVIIII we turn to systems with the nontrivial 
vacuum, using the Schrodinger equation. In Sec. IIXI we describe the path integral approach to the large- Ac models 
of the baryon wave function. In Sec. 0we check the agreement of the results based on the Schrodinger equation and 
on the path integral approach. 

II. MODELS 
A. Choice of models 

As was explained in the introduction, our aim is to check conjectures about the structure of the 1/Nc expansion 
for the baryon wave function which were put forward in Refs. 0|, 0|. We want to test these conjectures in simple 
large- Ac models. Our choice of these models is determined by the balance between the intention to preserve the main 
features of the 1 /Ac expansion in QCD and the possibility to solve the model. As was explained in the introduction, 
the aspects of the phenomenological relevance will play a secondary role (if any) in our choice of models. 

The minimal constraints on the models include two conditions: the model must have 

1) quarks degrees of freedom, 

2) SU{Nc) symmetry. 

The quark degrees of freedom are understood here in a rather weak sense. Neither full-fledged quantum field theory 
nor complete quantum mechanics is needed for our aims. For example, we can work with models ignoring the space 
motion of quarks where the dynamics of quarks is described only by color and some other discrete quantum numbers 
(e.g. spin and flavor). For the construction of the 1/Ac expansion we need the global SU{Nc) symmetry but there is 
no need in the local gauge invariance. 

Now let us turn to the dynamics of the models. Most of our work will be done in the Hamiltonian approach based 
on the analysis of the stationary Schrodinger equation. In principle, we can include both quarks and antiquarks in 
our models using the old Dirac picture. In this case the nontrivial physical vacuum will contain the "Dirac" sea of 
quarks above the vacuum and we must solve the Schrodinger equation both for the vacuum |0) and for the baryon 
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\B) 

i7|0) =fvac|0) , (2.1) 

H\B)=£b\B). (2.2) 
The dynamics of our models will be formulated in terms of the fermionic annihilation and creation operators a„cj 

where 

c= l,2,...,7Ve (2.4) 

is the color index and n — 1, . . . , K is some index which may have the meaning of spin, isospin, etc. 
We will be interested in the baryon wave function given by the transition matrix element 

{0\Y[a,njB). (2.5) 

where Qmc are quark annihilation operators. We can define the generating functional for this wave function by analogy 
with Eq. lfT3|) 

1>s(5) = (0|n fe.9'»«"^) 1^)- (2-6) 

c— 1 \ m / 

At large Nc we expect the exponential behavior (jl.6|l 

^B{g)'''^'^exp[N,W{g)]. (2.7) 

B. Hamiltonian 

We will work with models described by the Hamiltonian 

nin2n3"4 \c' = l / \c=l / nin2 \c=l / 

assuming that 

^l"2»3"4 ~ ^3n4nin2 • (2-9) 

We we will be interested in color singlet states. The color singlet states can contain MNc "quarks" where M is an 
integer number. For any fixed M we can study the problem of the lowest color-singlet state containing MNc fermions. 

In our models of the baryon function we associate notation AI with the vacuum whereas the baryon will correspond 
to M + 1 so that one has a nonzero matrix element l|2.5|l which will be interpreted as the baryon wave function. 

C. Hartree equation 

Below we will study the limit of large at fixed values of the parameters Vnin2nzni and L„j^„2 in the Hamiltonian 
(|2.8() . It is well known (see e.g. 0ill3) that in the leading order of the 1/Nc expansion the energy of the ground 
state can be found by solving the Hartree equation 
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where the single particle Hamiltonian /irii„2 given by 

M 

hnin2 = ^nin2n3ni4>ns,'t'ni + -^nina (2-11) 

and the single particle eigenstates are normalized by the condition 

^0;C=^"- (2.12) 

n 

In the leading order of the expansion the energy of the corresponding state is 

E = N,Eo + 0{N°) , (2.13) 

where Eq is given by 

^ M M 
-E'O = 2 X! "2 713714 '^"I'^ns'^ra'/'ni + X! X! ■^"l»2 ^nt 0n2 ' (2-14) 

r,S— 1 711712773714 S — 1 71i7l2 

In the case of Hamiltonians with Ln^n-i = the expression for Eq can be simplified using Eq. (|2.1U|) : 

1 

Eo = ^Y.^^ (ifi7n7i-0). (2.15) 

s=l 

Below we will see that the knowledge of the solutions of the Hartree equations for the vacuum |0) and for the baryon 
\B) is not sufficient for the calculation of the generating functional W{g) appearing in Eq. 12.7|l . One has to derive a 
special equation for W{g). 



D. Time-dependent Hartree equation 



The time-dependent Hartree equation (TDHE) appears in various problems of the traditional many-body physics 
(where its Hartree-Fock modification is usually discussed, see e.g. 0|). In the case of the Hamiltonian H2.8|) TDHE 
has the form 

'^CW^E'^^i'^^WCW, (2.16) 

712 

where 

1712 713 714 711712 . 

(2.17) 

r— 1 71304 

Although the context of the traditional problems of the many-body physics where TDHE appears differs from our 
problem of the calculation of the functional W{g), we will see in Sec. IVIHI that the functional W{g) allows for an 
interesting representation in terms of the solutions of the TDHE (|2.1()l) . (|2.17|) . The possibility of the Hamiltonian 
interpretation of the Hartree-Fock equations was discussed in Ref. ^.2_7.] . This Hamiltonian interpretation plays an 
important role in our analysis. 



III. MODELS WITH THE TRIVIAL VACUUM 



A. Simplest models 

As was explained in Sec. Ill Al we want to study the baryon wave function in models described by the Hamiltonian 
(|2.8|) using the states containing NcM and Nc{M + 1) fermions for the vacuum and baryon, respectively. 

We want to start from the simplest case M = (the generalization to arbitrary M will be considered in Sees. 
IVIIIllX)l . In these M = models the physical vacuum coincides with the bare one, and the baryon consists of Nc 
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quarks only. We will refer to the M = models as models with the trivial vacuum. The main advantage of these 
models is that the solution ^'bar of the Schrodinger equation 1)2.2(1 directly gives us the baryon wave function (|2.5I) . 

Our first aim is to define the generating functional $B(g) for the baryons in these models and to rewrite the 
stationary Schrodinger equation in terms of ^sig)- This can be easily done using the well known holomorphic 
representation (also known as boson representation in the context of the many-body physics ,.23j ) . 

In our simple models "baryons" are described by states containing Nc quarks: 

^ ]V~! ^ ^ ^ci...cjv<:''/'ni...njv,'^niV<3Cjv^ ■ ■ -""icJO) • (3-1) 

^' 71-1712 -..nj^^ CiC2...Cj^^ 

Here ipni...nK^ are coefficients. Due to the antisymmetry of fermionic operators coefficients ?/'ni...njv^ are completely 
symmetric in all indices ni . . . nj^^. 



B. Generating function for baryon wave function 



In our simple model the baryon wave function (|2.5|l coincides with the coefficients „ of the decomposition 

<...n„, =(0Kl...4„,e„JV'''). (3.2) 

Now we can construct the model analog of the generating functional (|1.5|) describing the baryon wave function: 

= E V'rfi...n„,5ni • ■ • 5««, = ^ (0|a+ ^ . ■ • a^„^c«J'/'^)5ni • ■ • • (3-3) 

nin2---njVc nin2...njVc 

Note that function ^sig) l|3.3(l is a homogeneous polynomial of degree Nc'. 

$b(A.9) = A^=$s(5). (3.4) 

The correspondence between the states \ip'^) and hmctions $_B(.g) is useful not only for the description of the 
eigenstates of the Hamiltonian (|2.8|) but for any color singlet states containing Nc quarks. Let us consider operator 



E ^mc^nc (3.5) 



c=l 

on the state In terms of the ^nig) representation for states 1-0^) 

IV'^) - ^B{g) (3.6) 

operator T,„„ takes the form 



C. Schrodinger equation for ^B{g) 



Using the g representation H3.7|l for the operators Tmn, we can rewrite the Hamiltonian (|2.8|l in the form 

^ = ^ E {gn. ^) (5«3 + E (5-. ^) . (3.8) 



ni n2n3n4 



Then the Schrodinger equation 

H\i;^) = SbH"") (3.9) 
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becomes 



2N, 



nin2n3n4 



5"i 



(3.10) 



IV. LARGE-A^c LIMIT 
A. Two approaches to the l/N^ expansion 

In this section we want to study the large- A'c hmit in the simplest models of baryons introduced in Sec. IIIII As 
was mentioned in the introduction, one can use two approaches to the large- iVc limit in fermion systems: 

1) approach based on the single-particle density matrix obeying the Hartree equation, 

2) approach based on the baryon wave function described in terms of the generating functional 

In principle these two methods are equivalent. For example, they lead to the same l/iV^ expansion for the spectrum 
of states. The second method based on the functional <I'_b((?) provides more information. However, one has to pay some 
price for the additional information: the large- iVc equations derived for the functional $5(17) are more complicated 
than the Hartree equation for the density matrix. 

Below we consider both approaches and demonstrate their equivalence. 



B. Large Nc Hmit for functionals ^sig) 

We want to study the limit of $B(.g) at large Nc- According to Eq. (|3.4|l we have 

W{Xg) = W{g) + \n\, 

so that 

T.Sn^^^Big) = 0. 



(4.1) 
(4.2) 



(4.3) 
(4.4) 



Let us insert Eq. H1.7|l into the Schrodinger equation (|3.1U|) . Taking into account the 1 /Nc expansion for the energy 



Sb = NcEo + AEb + 0{N-^) , 



we obtain in the leading order 

In the next order of the l/N^ expansion we find 

dW{g) 



dW{g) 
dgn2 



dWig) 



' dgn2 



9n3 ' 



dgn 



dlnAsig) 



dgn 



ogn2 



dWig) \ 
dgm ) 



5«r 



dW{g) 
dgn2 



Eo. 



3m ■ 



AEf 



dlnAsig) 



dg-n 



(4.5) 



(4.6) 



(4.7) 
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C. Hartree equation 

As was explained in Sec. Ill CI the energy of the lowest baryon state can be described in the leading order of the 
\/Nc expansion by the Hartree equations H2.10|l . (|2.11|) . Our case corresponds to Af = 1 in these equations: 



r!2 \n3n4 



3ni(l>n.^4>ni + -^niri2 j 4>n2 = £O0ni , (4.8) 

(4.9) 



where (/>„ is the single-particle wave function of the level occupied with quarks. 
The leading order contribution NcEq to the energy £b H4.5|) is given by Eq. H2.14|l 



in2 '/'ni </^n2 • 

(4.10) 



2 

nira2"3n4 nin2 

D. Agreement between the Hartree equation and the equation for W{ij>) 

Now we can compare equation (|4.6(l for W{g) with the expression (|4.10(l for the Hartree energy Eq. We see that 
the two expressions agree if we take 

, dW(g) 

K = , 0„ = g„ . (4.11) 

ogn 

In other words, if we know function W{g), then the Hartree solution (/)„ can be found by solving the equation 



9n = — =^ <Pn = 5„ ■ (4-12) 



Note that the normalization condition (|4.9|) is automatically satisfied 

W(., 
dgn 



EC^" = E5«^ = 1 (4-13) 



due to the property H4.3|l of the function W{g). 



E. Normalization integral 

Let us consider two states (|3.1|l 

|V'^'^>= E <^...n„.<,^.---<i|0) (fc = l,2) 



nin2...niv^ 



described by functions H3.3|l 



Then 



\r ir / / y iii2...iiVc «i 



^(2) 
ii2-..ijv„ 



Jll2---*JVc 



1 \ V 5 a 



1 d d 



N^^ J dgdg* exp {-N, 5.9*) [^^^^ (s)] * [^^'^ (5)] 



/d5d5*exp(-iV,^^5,5*) 



(4.14) 



(4.15) 



(4.16) 
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At large Nc we can use the asymptotic expression (|1.7|) for As a result, we find 

, (1), (2). e^-jV-^+-^ Jdgdg*exp{N4-g,g* + [Vt/(g)]* + Ty(g)} } [^(^(g)] * A(^)(g) 
V2^c Idgdg*e^p{-N,J:^9^9:) ■ ' 

This integral can be computed using the saddle point method. The saddle point equation is nothing else Eq. (|4.12|) 
with the solution g„ = 0* given by the Hartree equation H4.8|l . Thus we find with the exponential accuracy 

(V'WlV't^)) ~exp{-iVe{[VK(0*)]* ■ (4.18) 

We cannot use the saddle point method beyond this exponential accuracy because of the noncommutativity of the 
limits 



which is discussed in Sec. IV Bl 
Taking a normalizable state 



we conclude from (|4.18() that 



Nc^oo (4.19) 



^(1)^^(2)^^^ (4.20) 
(V'IV') = 1, (4.21) 



ReW^(0*)=O. (4.22) 

V. LARGE Nc LIMIT AND RPA EQUATIONS 
A. Beyond the leading order of the 1/Nc expansion 

In the previous section we have established a connection between the traditional description of the large- iVc systems 
in terms of the Hartree equation H4.8|l and our equation H4.6|l for the functional W{g). As is well known, the Hartree 
equation determines only the leading 0{Nc) part of the energy of the states. If one is interested in the 0{N^) 
corrections, then one has to solve the random phase approximation (RPA) equations 

On the other hand, in the problem of the baryon wave function the analysis of the next-leading-order of the 
1/A^c expansion is based on Eq. 14.7(1 for the functional As(g). Therefore one can expect that there must be some 
connection between the RPA equations and Eq. (|4.7|) . From our experience with the leading order equations [the 
Hartree equation H4.8() and Eq. 1(4.61) for VF(g)] we know that this connection may be rather nontrivial. Indeed, we 
will see that the RPA equation really appears in the analysis of the 1/A^c expansion for the basic functional $B(g) 
but in a modified version of the large- iVj, limit when the argument g of the functional $_B(g) changes with Nc in a 
special way. This new version of the large- iVc limit is discussed in Sec. IV Bl In Sec. IVDI we turn to the derivation of 
the RPA equations corresponding to the modified large- A^c limit. 

B. Two versions of the large-TVc limit 

Eq. 1(1.7(1 describes the asymptotic behavior of ^(g) in the large- A^c limit when g is fixed. 

g = const, Nc ^ oo . (5.1) 

But we can also study another limit 

A^c^oo, gn = (j)n + —=An , A„ = const . (5.2) 

when g„ approaches the Hartree solution in the limit Nc — > oo. In addition let us impose the condition on A„ 

^CA„ = 0. (5.3) 
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In the limit H5.2|l we have the asymptotic behavior 



$5 I 0„ + -i=A„) N^^Cb{A^) exp [N,W{c^)] (5.4) 



different from Eq. (|1.7|l . Taking A„ = and comparing the asymptotic expression (|5.4|) with Eq. p.7|) . we obtain 

C7b(O) = ^b(0). (5.5) 
Replacing A„ — > A^^A„ in Eq. (|5.4|l and using relation (|3.4(l . we find 

^£^=«iV,^BCB(A-iA„)exp{A^4W^(0)+lnA]} . (5.6) 

C. Special choice of basis 

It is convenient to choose the basis diagonalizing the Hartree Hamiltonian such that 

C - Sn • (5.7) 
In models with one occupied level we will label this level with the index s = and use the short notation 

^n^c^l^Sl. (5.8) 

Then the Hartree Hamiltonian H2.11|l becomes 

hmn — VnmOO + -^mn (5-9) 

and the Hartree equation (|4.8I) takes the form 

VmnOO + -trim = Sm^mn ■ (5.10) 

In particular, 

£o — ^0000 + -^00 • (5-11) 

The Hartree energy is 

Eq = -Vbooo + ioo • (5-12) 
We introduce a special notation for the components .g„ with n > 0: 

gn ffn , ffo = (5.13) 

so that 

9n = go<l>7i + gn , (5.14) 
*b(5) = *i3(5o,5) ■ (5.15) 

Now wc find from Eq. (|5.6(l 

l^B (<?)]^„=A„/vW: = *B {904". + ^) A^e^ooA^^t ^.uc^f^g^^A^) exp{iV, [W{^) + In go]} • (5.16) 
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Therefore 



d_ 



Nc — >oo,Ati —const 

dgo 



9o£r i^B (5)]^„=A„/v^ .v.^™-..... g^_N::-CB{9o'A,,) exp{7V, [W{c^) + Inyo]} 



{N^^ exp {N, [Wi^) + In go]}} 1 



d 



Cb(A„) 



(5.17) 



On the LHS the differential operator god /dgo acts on $b (5) whereas on the RHS we have the corresponding operator 
acting on Cb(A„). Equation (|5.17|) means that the following A representation is valid for the operator god/dgo- 



Similarly we find 



d 



go 



dgm 

gn^ > A„ 

ogm 



d 



9A„ 
d 



(5.18) 

(5.19) 

(5.20) 
(5.21) 



D. RPA equations from the large- A^c limit 

Now we can apply the results H5.18|l - H5.21|l to the Schrodinger equation (|3.10() . The terms of order Nc generated 
by god/dgo give the Hartree energy (|5.12|) . The terms of order \//Vc cancel because of the equation H5.1()|l . The terms 
of order 7V° lead to the equation 



E 

m,n>0 



Vo, 



A„ 



d 



dA„ 



1 ( ~ d 



:^A, 



1 ~ ~ ~ d 

2 dAn 



(VboOO + -^00) ^mn^m 



d 



9A„ 



Cs(A) 



AEBC'BiA). 



This equation can be simplified using Eq. H5.11|l : 

d 



E 



m,7i>0 



dA,, 



It/ fx ^ 

^KnOOn (^A„.^ . ^^^^ 



9 X 

r^Ar. 



2*'°™°"9A™9A„ 



1 ~ ~ ~ d 

2 dAr, 



£oS A 



d 



dAr. 



Cb{A) ^ AEbCb{A) . 



Let us introduce operators 



(5.22) 



(5.23) 



dAr. 



6+ — A 



(5.24) 
(5.25) 



Then 



E 

m,n>0 



[6m, 5+] = Sr, 



+ -;:VmOnob — SO^mnb^bn 



Cb = AEbCb 



(5.26) 



(5.27) 
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Note that the 6„ vacuum |0) 

6™|0)=0 (5.28) 
corresponds in the A representation l|5.24|l . (|5.25|l to the wave function 

*(A) = 1. (5.29) 

Thus we have the Hamihonian 



hb.pa = 



m,n>0 



jj_ I bjn \ I VornOn VnOQm + VoOnm + -^nm — ^oSmn 1 ( 



2 ''-^ \ bm J \ VmOOn + VoOmn + Lmn — Eo^mn VmQnO J \ b^i 



m>0 m>0 



Now we find 



where 



n>0 



Matrices a, /3 obey the conditions 



- ^ ^ [{Voomm + L.m,n) - Sq] . (5.30) 

m>0 

According to Eqs. H2.9|l and H5.10|l we have 

^OOnm ^^ ^nm ^ ^m^nm i (5.31) 

^ ^ (Vbomm + Lmva) = ^ ^ £?n . (5.32) 



^HPA - ^ E ( )'^^ ( ) - ^ E - -o) , (5.33) 

m,n>0 ^ ^ ^ ^ m>0 

■J^ — f K)mOn KiOOm + (^m — £o) <5mn j j-g g^-j 

This is nothing else but the large- A^c version of the RPA Haniiltonian. 

E. Diagonalization of the RPA Hamiltonian 

We can diagonalize the RPA Hamiltonian (|5.33|) using the Bogolyubov transformation 

UbmU^^ = E {^rnnbn + Pmnb^) , (5-35) 
n>0 

Ub+U~' - J2 i^'^nbn + a*„,X) , (5.36) 

n>0 



U-'b^U - J2 «r,^bn - Pnmbt ) , (5.37) 
U-^JJ = E i-f^nrnbn + «„«&„) ■ (5.38) 



aa+ - P(3+ = 1 , (5.39) 
aP^ = (3a^ , (5.40) 
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a+a - P^f3* = 1 
Eqs. H5.;-{9|l and (|5.4()|l can be rewritten in the form 



Introducing the notation 



we can write 



f a 




( a+ 








1-/3+ 





a P 
(3* a* 



u\ II If/- 



s 



b 

b+ 



Now we can diagonalize the RPA Hamiltonian H5.33|l using Eq. H5.45|l 



- E 

m,n>0 



Let us choose S so that 



where is a diagonal matrix 



^ m>0 



n 



Q = diag(r2„) . 



Then 



UHrpaU ^ = J E (^™^™ + ^™^") " ('^"^ ~ ^o)] ■ 



m>0 



(5.41) 
(5.42) 



(5.43) 



(5.44) 



(5.45) 



(5.46) 



(5.47) 



(5.48) 



(5.49) 



The spectrum of this Hamiltonian gives the 0{N^) contribution AEb to the l/Nc expansion of the total energy (|4.5|l : 



m>0 



The corresponding eigenstates of i?RPA are according to Eqs. H5.38|l and (|5.49(l 



n(^ 



k ) 



m>0 



C/-i|0) 



Here |0) is the bm vacuum H5.28|l and the ground state of -//rpa is 

IOrpa) =C/"'|0). 

Using Eqs. ifCT]! and 1^231, we find 

= U-^b„M - {U-^bmU) U~'\0) = Krnbn " Pnmbt) C/"'|0) 



n>0 



Combining Eqs. 15.52|1 and (|5.53|) . we obtain 

E i^nmbn - /3nm&+) |Orpa) = . 



(5.50) 



(5.51) 



(5.52) 



(5.53) 



(5.54) 



n>0 
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F. Calculation of functions Cb(A) 

Functions Cb(A) were defined via the asymptotic expression (|5.4|1 . Now let us compute these functions. We start 
from the function Co(A) corresponding to the ground state |Orpa)- The ground state |Orpa) obeys Eq. (|5.54|l . In 
the A representation H5.24|l . (|5.25|l . Eq. H5.54|l takes the form 



E 

ri>0 



d 



/3„™A„ Co(A) = 0. 



The solution of this equation is 



where 



Co(A) = AAexp 



i(AXA) 



(AX A) = ^ KaXmnK. 



(5.55) 



(5.56) 



(5.57) 



m , Ti > 



and A/" is a normaUzation factor. 
According to Eq. H5.42|l we have 



/3 (a*)-^ = V^. 



This means that matrix 

J van 

is symmetric 

The excited states H5.51|l are described in the A representation H5.24|l . ()5.25f) by the wave functions 



d 



Cb{A) = ( ^ -^*mk-^ + amk'^m ] exp 



k \m>0 

Here Pb{A) is a polynomial of degree 



At large A„ we have according to Eq. 



i(AXA) 



7VPB(A)exp 



l(AXA) 



k 



where 



PB(A)^^n[(-'3^^^ 



{-P+XK)^^-Y.PrakXmnK. 



Using Eqs. H5.60|l and H5.41|l . we find 

= P+p {a*y^ = (a^a* - l) (a*)"' = - (a*)"' 



(5.58) 



(5.59) 



(5.60) 



(5.61) 



(5.62) 



(5.63) 



(5.64) 



(5.65) 



(5.66) 
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G. Matching two types of the 1/Nc expansion 



In Sec. IV Bl we have described two different asymptotic expressions for the large- A^c hmit. These two expressions 
are vaUd in different regions. However, we can match the two expressions in their common overlap area. 
Let us insert Eq. (|5.t)2|) into Eq. (|5.4|l 



A 

TWr 



= JVPB{A)exp 



exp [N,W{^)] 



This expression must agree with with Eq. H1.7|l 



exp 



NcW <j> + 



A 



, A„ = oivK) 



(5.67) 



(5.68) 



in the region 

l<A<v/^- (5.69) 
In this region only the leading term of degree (|5.64|) survives in the polynomial Pb(A) so that Eq. H5.67|l becomes 



A 

Wr. 



exp 



i(AXA) 



exp [NrW{^)] 



Now we want to check the general factorization relation H1.18|) 

Ab (g) = AW(5)n [^^(5)]"'" = A^°\9)l[[Ak{~g) 

k k 

If this factorization holds, then Eq. H5.68|l takes the form 

'a 



X exp NcW {(j) - 
Due to Eq. H4.11|l and condition H5.3|l . 



A, 



k L 

A„ - O(yiV^) 



= 0, 



the variation W ~\- S(fj has no linear term: 

w (0- 

Thus 



$B + 



A 



i(°)(o)n 



Ak 



exp 



Inserting this into Eq. H5.70|l . we see that 
= 7V,''«i(o)(0)n 



Ak 



^ exp 



i(AXA) 



exp [NcW{cj))] 



exp 



(5.70) 



(5.71) 



(5.72) 



(5.73) 



(5.74) 



(5.75) 



(5.76) 
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Now we see that 

i(")(0)=AA7Vc"'^"+^ (5.77) 
Ak (0 + (50) = (-I3+X 5$) ^ , (5.78) 

W'^^\54),5^) = {54>X5$) . (5.79) 

H. Asymptotic behavior in the vicinity of the Hartree solution 

Inserting Eq. (|5.79|) into Eq. H5.74|l . we find 

W [<j) + (5^) = VF(0) + i [54: X5!j^+0 ((50^) . (5.80) 

Combining this result with Eq. H4.1|l . we obtain 

W {\(t) + 5^^ = W [<j) + X-^5^^ +lnA W{<t>)+\nX + ^ [54) X 5$^ + O {s^^^ . (5.81) 

Taking into account Eq. (|5.7|) . we conclude that if 

l5m|«.9o (m>0) (5.82) 

then 



^ V gmXmngn + 0{ (^-^\^\ . (5.83) 

?o) ™^„>0 ^0 ^ / 



W^(ff) = l^(,^)+ln5o+ , , 
2 (50) 

In Sec. IIV Bl we derived Eq. (|4.t)|) for the function W (g). This is a partial differential equation which has many 
solutions. The asymptotic expansion (|5.83|) plays the role of the boundary condition for this differential equation. 
The matrix X^n appearing in this boundary condition is given by Eq. I|5.60|) . 

In the same way, using Eq. (|4.2(l . we can rewrite relation (|5.78|l in the form 



Ak{g)=-- P*mkXrangn + 0{(^ 



2\ 

.9-m>0 



(5.84) 



go 

This asymptotic behavior of Ak [g) should be used as the boundary condition for the differential equation 1)4. 7|l . 

VI. CLASSICAL DYNAMICS OF LARGE-Af^ SYSTEMS 
A. Hamilton— Jacobi equation 

According to our previous results, function W{g) can be computed by solving the differential equation ()4.f)|l with 
the boundary condition (|5.83|) . This representation for W{g) is not quite convenient since one has to deal with the 
partial differential equation H4.6|l . In this section we show how the problem can be reduced to the analysis of ordinary 
differential equations. This reduction can be done using the well-known Hamilton-Jacobi method. 

Eq. H4.6|l has the structure of the classical Hamilton-Jacobi equation. In order to make this structure explicit, we 
change the notation 

g^q. (6.1) 
W{g) S{q) . (6.2) 
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Then Eq. H4.6|l takes the form 
where the Hamiltonian is 

in2"3n4 (gniPna) (g^aPru) + X! -^"i"^ ilniPni) ■ (6-4) 

One should keep in mind that g can be complex and W{g) is an analytical function of g (which may have singular 
and branch points). Therefore our classical mechanics is complex. 

As is well known, the solutions of the Hamilton-Jacobi equation can be constructed in terms of the action for the 
trajectories obeying the Hamiltonian equations 

dpn _ dH{p,q) dqn _ dH{p, q) 

dt dqn ' dt dpn 

Solution S{q) of the Hamilton-Jacobi equation H6.3|l describes the set of trajectories covering the coordinate space. 
All these trajectories have the same energy Eq and obey the condition 

Pn = ■ (6.6) 

oqn 



B. static and time-dependent Hartree equations 

If we change the notation 

071 ^ Pn , (t>n^ Qn (6.7) 

in the Hartree equation 14.8|l and in the normalization condition H4.9(l then we obtain 

d 

-7^H{p,q) ^ eoqn, (6.8) 
opn 

d 

j^H{p, q) = eopn , (6.9) 
oqn 

^PnQn^^, (6.10) 

n 

Pn = qn ■ (6.11) 

Let us stress that we consider p„ and qn as independent complex variables of our phase space. Conditions (|6.8|) - (|6.11|) 
define a point in this space. At this point the variables p„ and qn are complex conjugate. But generally speaking 
Pn ^ Qn- 

Strictly speaking, Eqs. I|().8|l - I|t).ll|l have more than one solution. For example, these equations are invariant under 
the phase transformation 

Pn e^'^pn , ^71 e'^^qn ■ (6.12) 

In principle, there can be additional symmetries which increase the number of solutions. 
If one chooses the basis (|5.7f) . then the solution of the Hartree equations (|6.8|l - Ht).ll|l is 

Pti = qn = SnO ■ (6.13) 

Parameter Eq is given by Eq. H5.11|l . 

The change of notation 16.7|l in the Hamilton equations 16.5(1 transforms them into the form 

= - X! ^^nn(t)n , = ^ K^hmn , (6.14) 

71 m 

where 

hmn — y ] Vmnijff'i'f'j + -^7Ti7i ■ (6.15) 

These equations can be interpreted as a Euclidean version of TDHE (|2.1(j|) . (|2.17|) . It should be stressed that the 
variables (j)n and (/>* are not complex conjugate in equations (|6.14|) in contrast to the standard TDHE. 
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C. Asymptotic conditions 

In the vicinity of the "Hartree point" H6.8|l - H6.11|l . the Hamilton equations for trajectories can be approximated by 

dqn ^^^p^ ^ ^^^^ ^ {e.l6) 



dt dpn ' 

dpn 9 

dt dqn ' 



- .J_H{p,q)^-S0Pn (6.17) 



with the asymptotic solutions at t — > — oo 



qn{t) SnQ exp [eo [t - t)] , (6.18) 

Pn(i)*^=°°<5noheo(t-T)] . (6.19) 
According to Eq. (|6.4|l this asymptotic behavior corresponds to the Hartree energy Eq (|5.12|I 

H{p,q)^Eo. (6.20) 

D. Integral of motion 'Y^^PnQn 

The Hamiltonian H6.4|l and equations of motion (|6.5|l are invariant under the transformations 

Pn A" Vri , qn ^ Mn ■ (6.21) 

Therefore 

Pn^ g„— )i/ = 0. (6.22) 

OPn Oqn J 

This identity can be rewritten in terms of the Poisson bracket 



J2Pnqn,H\=0. (6.23) 

. n ) 



Thus we have the integral of motion 

'^Pn.qn = const . (6.24) 

n 

For the trajectories with the asymptotic behavior H6.18|l . (|6.19|) we have 

^P«'7n = l. (6.25) 

n 

E. Boundary condition 

Using notation 1)6. l|l . ()6.2|l . we can rewrite the asymptotic expansion H5.83|l in the form 



S{q)= Siq")\^,^^^^^+lnqo + l ^ X^^^^^ + O 



Qni Qn 



2 qo qo 

m,n>0 



qm>o 



qo 



(6.26) 



This expression is valid for 

\qm\<^qo (™>0). (6.27) 
Inserting this decomposition into Eq. (|6.6|l . we find 



Puiq) = = -Sno if \qm\ « 90 (m > 0) . (6.28) 

oqn qo 
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This agrees with the asymptotic behavior of trajectories H6.18|l . (|6.19(l a,t t ^ —oo. 

Now we understand that the action S{q) is associated with the configuration of trajectories p{t), q{t) which start 
at t ^ — oo at points p, q obeying Eqs. H6.8|l - H6.1()(l with the asymptotic behavior Hfi.1811 . Hfi.l9|l . All these trajectories 
have the same energy Eq and the same integral of motion It). 2511 . 

Let p{t), q[t) be one of these trajectories with energy Eq: 



H[p{t),q{t)]^Eo. 

Let us assume that at t = this trajectory passes through the point p'-*'^ , (7^"-' : 

p(0)=p("), g(0) = g("), 



(6.29) 



(6.30) 



and at i ^ —00 it has the asymptotic behavior (|6.18|) . H6.19|) with some 

r-T(g(")) 

depending on q^: 

gn(<) <5„oexp jso t-T{q^"^) I , 

Note that at i ^ —00 we have along this trajectory 

L = ^Pnqn - H{p, q) eo - ^^o • 



This means that the integral 



f dt{L + Eo 

J —00 



is convergent at t ^ — oo. Now let us define the function 



S{q^"^)= / di(i + i?o-e„)-sor(gW) 



and show that that this function obeys the differential equation (|6.3|l . 
According to Eq. H6.29|l we have 

L + Eq ^ L + H ^y^^Pnq„ 



Therefore 



dt 



^P,i{t)qn{t) - £0 



Hence 



/■O 

5(g(°))= hm / dty^pn{t)qn{t) + SoT 

T^-°°JT+Tiqm) „ 



= lim 

T— >-oo 



/ y^Pndqn + epT . 

J5„oexp(£oT) „ 



(6.31) 

(6.32) 

(6.33) 
(6.34) 

(6.35) 

(6.36) 
(6.37) 



(6.38) 



(6.39) 



(6.40) 
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Taking some large but fixed |T| on the RHS, we see that this definition of S'(g(''^) has the standard property 

(6 41) 

where pn"^ is the momentum Pn{t) taken at t = according to Eq. (|6.30() . Inserting Eq. H6.41|l into Eq. (|6.31|) and 
taking into account that the point g'^"-' is arbitrary, we find that the function S{q) defined by Eq. H6.37(l obeys the 
Hamihon-Jacobi equation ^<6.'6\ . 

Since the definition Ht).37|l of S{q) is based on the trajectories with the asymptotic behavior H6.33|l . (|6.34|l . we 
automatically have the property H6.28|l . Therefore function S{q) defined by Eq. (|6.37|) also obeys the boundary 
condition (|6.26(l . Since both differential equation l|6.3|l and the boundary condition H6.26|l are satisfied, we conclude 
that Eq. H6.37(l gives a correct representation for the functional 

(6.42) 

which determines the large- A^c behavior (|1.7|l . 

The construction described in this section leads us to the following algorithm of the calculation of W{g). 

1) Solve equations of motion H6.5|l imposing boundary condition 

g(0) = g (6.43) 
with given g ()6.60|1 and boundary conditions H6.33(l , l|6.34|l 

qn{t) *^=°° <5„o exp {eo [t ~ T{g)]) , (6.44) 
Pn{t) <5„o {-£0 [t - Tig)]} . (6.45) 

with unknown T{g). In this way one finds the trajectory and T{g). 

2) Now W{g) is given by Eqs. 



W{g) = / dt 



— OO 



^Pn{t)qn{t) 



£0 



eor(5) . (6.46) 



F. Case Lmn = 

Let us consider the special case when 

Lran = (6.47) 

in the Hamiltonian H6.4|) . In this case the Hartree equations (|5.10p . H5.12|l lead to the simple relation 

- ^£0 ■ (6.48) 



In the case ()6.47f) the Hamiltonian 1)6. 4|l is quadratic in p so that 

H{p, 

dpn 

Combining this with Eq. H6.29|l . we find 



Y.Vr.'-^ = 2Hip,q). (6.49) 



5lM=IIP"7^ = 2F(p,g) = 2£;o. (6.50) 

71 n ^ 

Taking into account Eq. (j6.48|l . we arrive at 

^Pn^n - £o = 2£;o - £o = 0. (6.51) 



25 



Now we insert this into Eq. H6.39|l 

5(g(°)) = -eor(<7(o)). (6.52) 
This ahows us to rewrite Eqs. H6.33(l and H6.34(l in the form 

quit) <5„o exp [eot + Siq^^^)] , (6.53) 

Pnit) Sno \-eot - 5((z("))l . (6.54) 



Thus in the case (j6.47|l the action can be directly read from the asymptotic behavior of qn{t), Pn{t) at t — > —oo. 
Remember that these trajectories are fixed by extra conditions (I6.30|l 

q{0) = (6.55) 

and KM . llOa 

=i?o = i£o. (6.56) 

Now we conclude from Eqs. 16.42|) . H6.53|l . (|6.54|) that in the case L,„„ = the functional W{g) is given by the 
equation 

W{g)^lnl{g) (ifL„,„ = 0), (6.57) 
where I{g) is the parameter of the trajectories obeying the asymptotic conditions 

qn{t)'^^°" SnoIig)e'"\ (6.58) 
p„(i)*-=-^e-°*, (6.59) 

9„(0)=.gn. (6.60) 

Let us summarize. In the case L„in = the calculation of the functional W{g) reduces to the following steps: 

1) Solve equations of motion 1)6. 5|l imposing boundary conditions H6.58|l - (|6.6Q(I with given g and with unknown 
I{g). In this way one finds I{g). 

2) Express W{g) via I{g) according to Eq. (|6.57(l . 

VII. EXAMPLE: ASYMMETRIC ROTATOR AS A LARGE- iVc SYSTEM 

A. Model 

In the above sections we have considered the general theory of fermionic systems described by operators a„c7 ^tc 
(|2.3|) with c = 1,2, Nc- Now we want to consider the simplest version of these models with only two values for 
the index n of a„c- It is convenient to interpret this two-valued index as the projection of spin 1/2. In other words, 
we want to consider the model described by the Hamiltonian 



H=^J2^-bJaJb+J2YaJa, (7.1) 
" ab a 

where 

We ^ 

T.o4cin),,,a,,. (7.2) 



Obviously this Hamiltonian belongs to the class 1)2. 8|l with parameters 

^in2n3n4 = ^ '^^^ ^ab i'''a)m no ('''b)n3n4 ' (7-3) 
ab 

Lnin2 ~ 77 ^ ^ (''"1)711 719 ■ (7-4) 
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On the other hand, Hamiltonian H7.1() is nothing else but the rotator described by the angular momentum Ja 

[Ja, Jb] = iSabcJc ■ (7.5) 

We will consider the case of the asymmetric rotator {lab 7^ I^ab) but the quantum number is still conserved: 

[J',H]^0 (7.6) 

and has the standard eigenvalues 

J2 = J(J+1). (7.7) 

We want to study the states made of Nc quarks and antisymmetric in color. The Fermi statistics and the color 
antisymmetry lead to the complete symmetry of the spin wave function so that we deal with 

J = ^ . (7.8) 

Thus the algebraic formulation of the problem of the baryon is obvious: take the sector with J — Nc/2 and diagonalize 
the Hamiltonian (|7.1|l in this sector. 



B. Direct semiclassical approach 

In principle, we can find the spectrum of low-lying states of the Hamiltonian l|7.1|l at large Nc solving the general 
Hartree equation (|4.8|) and diagonalizing the RPA Hamiltonian H5.33|l . However, the same results can be obtained by 
directly applying the semiclassical approximation to the large momentum J (|7.8|1 in the Hamiltonian (|7.1() . 

Instead of solving the Hartree equation (|4.8(l we can simply minimize the Hamiltonian H7.1(l considering it as a 
function of the classical momentum J. This minimization should be performed assuming the classical version of the 
constraint (|7.8|) . Introducing the Lagrange multiplier A, we arrive at the extremum problem for the ground state 
energy H4.5() in the leading order of the l/Nc expansion: 



Nr-En = min 
J 



^ {Kab - X5ab) JaJb + Yj ^""^-^ 



ab 



(7.9) 



This leads to the equation 



XSab) Jb + Ya^O. 



(7.10) 



If if — A is not degenerate, then 



Ja^-NcJ2[iK-Xr']abYb 



(7.11) 



Inserting this into Eq. I|7.8|l . we find 



{Y,{K-X)-^Y) = ] 



(7.12) 



This equation determines A. Knowing A, we can find J from Eq. H7.11|l . 

The leading order of the l/Nc expansion 1)4. 5|l for the energy of the lowest states is determined by the parameter 



1 



2N, 



■KabJaJb + YaJa 



-Y{K - X)-^K{K - X)-^Y - Y{K - Xy'^Y 



2 (K- A)2 



(7.13) 



In principle, we may also have other solutions of Eq. I|7.10(l corresponding to 



det{K- X) = 0. 



(7.14) 
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C. Special case 

Below we concentrate on the case when matrix Kab is diagonal 

Kab = KaSab (7.15) 

and 

^1=^2 = 0. (7.16) 

We also assume that 

>3>0, (7.17) 

ifi-X3 + 2y3>0, (7.18) 

+ (7.19) 

In this case Eq. H7.12|l yields 

- A = ±21-3 (7.20) 

and Eq. (|7.13|) results in 

So = (i^3 - 2A) = [K^ - 2 (if3 T 2F3)] . (7.21) 
According to Eq. (|7.17() the minimal energy corresponds to the upper sign. Thus 

We find from Eq. ifTTHl 



Eo^l (Ks - 4^3) . (7.22) 



Ji = J2=0, ^3 = -^. (7.23) 

One can show that under conditions H7.17|I " H7.19|I this solution gives the true semiclassical ground state and solutions 
of Eq. H7.14|l can be ignored. 

D. Hartree equation 



The semiclassical spin Ja has the following interpretation in terms of the solution of the Hartree equation 14.8|l 

2 



Ja^^Y.'l>*'ni^-)mn^rn. (7.24) 



mn 



Therefore the classical solution H7.23|l corresponds to 

(l)m = S,n2 = (^^^ ■ (7-25) 

In the Hartree equation (|4.8|l we have the single-particle Hamiltonian 

hmn2 = ^"l"2n3n4C3'^"4 + ^nin2 = ^ (--^^3 + 2^3) (t3)„i„2 " C^'^^) 

Hence the single-particle energy £o is 

eo = \ {K3 - 2Y3) . (7.27) 

It is easy to see that the general expression for the Hartree energy Eq 14.10|I reproduces the above result H7.22|l . The 
single nonoccupied Hartree state is 



(1) - 



(7.28) 



and its energy is 



ei = -eo = -^{K3-2Y3) . (7.29) 
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E. RPA equation 

In the general RPA matrix TZ H5.34|l . the indices m,n correspond to the nonoccupied eigenstates of the Hartree 
equation. In our case we have only one nonoccupied state H7.28|l . We label the occupied Hartree state with m = 
whereas the nonoccupied state is labeled by m = 1. As a result, matrix TZ (|5.34|) takes the form 



^1001 + £1 — £0 VioiO 
Using Eq. ()7.3|l. we find for the diagonal tensor Kat H7.15|l 

V^noi = V^ioio - J (^1 - K2) , (7.31) 
^1001 K2) . (7.32) 

With these coefficients Kiin2n3n4 and with expression (|7.29l) for Si — Eq, we obtain 

_}_{ K1-K2 K1 + K2 + 2 i-Ks + 2X3) \ .7 oo^ 

4{Ki + K2 + 2{~K3 + 2Y3) K1-K2 )■ ^' 

This leads to the Hamiltonian 15.33|l 

Hkpa = I [{Ki - K2) {aa + a+a+) + [Ki + K2 + 2 (-X3 + 2^3)] {a+a + aa+)] - ^ (-ifg + 2^3) • (7.34) 
o 4 

Performing Bogolyubov transformation 

«l/2 + ^-l/2 ^1/2 _ ^-1/2 

a= 6+ b+ (7.35) 

with 



K2-K3 + 2Y3 



K1-K3 + 2Y3 
we find 



(7.36) 



i^RPA = I {b+b + bb+) + i (Ks - 21^3) , (7.37) 

where 

w = i ViKi -K3 + 2Y3) {K2 -K3 + 2Y3) . (7.38) 
According to Eqs. ifT?^ and (|73S|) we have 

Ki= K3-2Y3 + 2k-^lo, (7.39) 

K2 = K3 - 2Y3 + 2kuj . (7.40) 

Comparing Eq. (|5.49|) with Eq. (|7.35|l . we see that 



l/.V2 + ,-l/2 ,l/2_,-l/2X 

" 2 - + y ^ ^ ^ 

and the only RPA excitation energy is 

Vti=u. (7.42) 
Comparing our result (|7.41|) for S with Eq. H5.45|l . we obtain 

ail = ^ ^-1/2 ^ (-7 43) 

= ^1/2 _ ^-1/2 . (7.44) 

Using Eq. (|5.t)U|) , we find 
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F. Method of generating functionals 

In the case of the Hamiltonian (|7.1|) . Eq. (|4.6() for the functional W{g) takes the form 

l^abJaJb + YaJa = Eq , 
ab a 

where 

Ja=9. (1)^ 



'ra\ dWjg) 



(7.46) 



(7.47) 



According to Eq. (|4.1|) function (31,52) has a nontrivial dependence only on the ration (71/(72 so that we can 
represent M^( (71,(72) in the form 



m5i,52) = /(^|) +^ln(5i52) 



(7.48) 



with some function /. It is easy to see that Eq. (|7.46|l leads to an ordinary differential equation for the function /. 
In order to solve this equation it is convenient to change the variables: 



In terms of these variables representation H7.48|l becomes 



W{gi,g2)^\R{0 + \w. 
A straightforward calculation of expressions H7.47|l shows that 



Jl = 


1 


(1 


J2 = 


i 


(1 




i.ai?(C) 

2^ ac 



' dC 

dm 

' dC 



In the case (|7?T5|) . (fTTCIl we find from Eq. (|7^ 



1 ^ 

-Y^aiJaf+Y^Js^Eo. 



Inserting Eqs. (|7.51(I - H7.53(I . we obtain 



2C 



(i + C)-(C-i)C 



dRiO 



1 2 



dC 

dRiO 



K2 

2C 



(i-C) + (C + i)C^ 

dC 



dC 



8r3C«^16i.o. 
dC 



(7.49) 

(7.50) 

(7.51) 
(7.52) 
(7.53) 

(7.54) 



(7.55) 



At the point g„ 



G. Agreement with the Hartree equation 
(|4.11|l corresponding to the solution of the Hartree equation we must have according to Eq. 



d 
dgn 



1, (I) UIm...) 



(7.56) 
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Taking into account Eq. (|7.25() . we see that this imposes the following constraint on the behavior of R{C) at small ( 

i?(C) const - In C + . . . (7.57) 

Let us find the next term of this small C expansion. Applying Eq. H5.83|l to our case and taking Xn from Eq. 
(|7.45|l . we obtain 

W{g) = +lng, + ^^^^-^ ( ilV + Q ( ( ^-^Y] . (7.58) 



According to Eq. (|4.22|) we have for normalizable states (with the appropriately chosen phase) 

W{(p)=0. (7.59) 

Therefore 

W (.) = In,. + l^Z^ ('-^X+o((iiy]. (7.60) 



2/tl/2 + ^-l/2 ^g^J ^^g^ 

Combining this with Eq. (|7.50() . we find 

m)'=-ir^c+2^^^^^c+o{e) . (7.61) 

H. Calculation of W{g) 

Eq. H7.55|l is a quadratic equation with respect to dR/dz which can be easily solved. Using expression (|7.22|) for 
Eq , we find from Eq. (|7.55|) 



dR ^ - [ri{l - e) + CC] ± 2C^ie ~ V') - yCC 

^dC ,7(l + C^)-(2?-C)C ■ ^ ' ' 
Here we have introduced the compact notation 

C = 4^3 , (7.63) 

e = + k)u; = ^{Ki+ K2) -K3 + 2Y3 , (7.64) 

r, ^ {k-' - k) CO ^ ^{Ki ~ K2) . (7.65) 
In order to fix the sign uncertainty on the RHS of H7.62|l . let us consider the limit C ~* 

^dR c^o , 2C 



-1 + ^ ± ve^) = -1 + 2c ^ ^ , (7.66) 

On the other hand, we find from Eq. (|7.61|) 

Comparing these two expressions, we conclude that we must have 

^1/2 _ ^^1/2 ^ _(^-i+k)±2 



This identity holds if we choose the upper sign in the numerator of the RHS. Now Eq. (|7.62|) takes the form 



(7.68) 
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One can easily integrate this equation. The integration constant is fixed by Eq. H7.61|l . The calculation of the integral 
yields 

R{C)=-\nC + F{C)-F{0), (7.70) 

where 

^ ^ vi In [u(C) - vi] - V2 In [ujQ - V2\ 



"^•^ ^ W^2\ ' ^'-"'^ 



"^^^ - \K,-K2\ • ^^-^"^^ 

Now we insert Eqs. (fTTUI) and ifT^ into Eq. lf73n|l 

1^(51,52) = In .92 + If |^(^|^ ^ - \f{Q). (7.74) 

Thus we have computed function W{gi,g2) which determines the exponential part of the large- A^c behavior l|1.7|) 
for the system with the Hamiltonian H7.1|l . One should keep in mind that function F{C,) has branch points. At small 
C this function is regular. According to Eq. H7.61|l 

F(C) F(0) + \!J2~S^U + O ie) . (7.75) 

Therefore for small I51/52I our result for W{gi,g2) is unambiguous. However, at larger \gi/g2\ one can meet cuts. In 
this case the determination of the relevant branch of W{gi, 32) will require a special analysis. 



VIII. MODELS WITH THE NONTRIVIAL VACUUM 



A. Color singlet states with MNc quarks 

Our previous work was devoted to models of the baryon wave function based on systems with the trivial vacuum. 
In such systems the baryon wave function comes directly from the solution of the Schrodinger equation for the baryon 
(|2.2(l . In other words, we worked with models which have quarks but no antiquarks. Within this class of models 
we could check the general properties of the large- iVc limit like the exponential large-iVc behavior l|1.7|l and the 
factorization of the preexponential factors Ab (jl.l8|l . 

Now we want to show that these general large-iVc properties also hold in more complicated models containing 
antiquark degrees of freedom in addition to quarks. It is convenient to formulate these models in terms of the old 
Dirac picture when both vacuum and baryons represented as a result of the occupation of the bare vacuum with 
quarks. In these models the physical vacuum contains MN^ quarks whereas the baryon is made of (M -I- 1)-/Vc quarks 
where M is some integer number. This allows us to define the quark wave function of baryons via the transition 
matrix element (|2.5ll . 

We want to start from the case M = 1 in order to concentrate on the idea of the method and to avoid unnecessary 
complications. The generalization for arbitrary M is discussed in Sec. IVIII ,71 in terms of the Schrodinger approach. 
In Sees. IIXI 1x1 the case of arbitrary M will be analyzed using the path integral method. It should be stressed that 
the case M = \ from which we want to start is rather nonphysical: the vacuum made of Nc quarks is a fermion for 
odd Nc whereas the baryon made of 2A''c quarks is always a boson. Nevertheless the M — 1 case is a good starting 
point for understanding the dynamics of models including quark-antiquark pairs. 
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B. Color singlet states with 2Nc quarks 



The color singlet states with 2Nc quarks above the hare vacuum |f2) can be represented in the form 

ikjk c=l c' = l 

Without any loss of generality we can assume that ^jjjj is symmetric under permutations of upper indices 
iii2 ■ . - iNa well as under permutations of lower indices jij2 . . . jN^'- 

„ili2...ijv^ _ p{ili2---iNc} (Q r)\ 



The Fermi statistics of quarks 



allows us to antisymmetrize B^j^^j^ '^jl^ with respect to ii ^ ji. 



aic4,c' = -atc'<c (8-3) 



2 y^hh-jN, - B^lJ2■■■3N,) ■ (8.4] 



We can continue this process antisymmetrizing in all pairs ik ^ jk- After that the original symmetry in {iii2 . . . ijVc} 
and in {jij2 . . • jN^} is lost but we still have the symmetry with respect to the permutations of any pair {imjm} with 
another pair {injn}- 

Note that our construction of tensors S!^!^ based on 

1) symmetrization in {11^2 . . . iN^i 

2) symmetrization in {jij2 • . -Jn^}, 

3) antisymmetrizations [iiji], [i2j2]v, [iNjN,] 

coincides with the construction of tensors associated with the rectangular Young tableau containing two rows and Nc 
columns. 

Assuming that the tensor Bl^^^'"^'^" obeys the above symmetry properties, let us construct the "generating func- 
tional" 

(7) = ^'hh'Xl ^nil 1^202 ■ ■ ■ IinJn, , (8-5) 

ikjk 

depending on the antisymmetric matrix "source" 

Irj = -In ■ (8.6) 

This construction generalizes the case of states with Nc quarks H3.3|l . We use the tilded notation $ in order to 
distinguish the functional $3(7) depending on the antisymmetric tensors 7ij from the functional $_B(.g) depending 
on "vectors" gi. According to Eqs. (|8.1|) . (|8.5() the functional ^b{i) corresponds to the transition matrix element 
between the baryon \B) and the hare vacuum |f2). On the other hand, the functional $s(g) is associated with the 
transition matrix element (|2.6|) between the baryon \B) and the physical vacuum |0). 

Now we have the correspondence between the states \B) (|8.1|l and functions $3(7) IjS.Sfl : 

\B) ^ $s(7) (8.7) 
Within this correspondence operators T^n H3.5() are mapped to 

d 



According to Eq. (|8.6(l tensors 7„j are antisymmetric. Therefore we must be careful about the normalization of the 
derivative d/djnj- Our choice is 



d 



pq 
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Obviously we have 

E^'J = 27Ve$B(7) . (8.10) 

Now we turn to the Schrodinger equation. For simphcity we will assume that 

L„„, = (8.11) 

in the Hamiltonian (|2.8|l . Using the 7 representation H8.8|l for this Hamiltonian, we can rewrite the Schrodinger 
equation 

H^b{i)^Sb^b{i) (8.12) 

in the form 

d \ d 



mnpq \ j 



f> ^ \Y."^pkl^]^B{l)=£B^B{l). (8.13) 



C. Large-TVc limit 



At large iVc we use the standard ansatz for the generating function corresponding to the 2N(, quark state: 

^5(7) = K^Ab{i) exp [iVeW^bar(7)] ■ (8-14) 

Inserting this ansatz into Eq. (|8.10|l . we find 

E^^^^^bar(7) = 2. (8.15) 

Combining the large- iVc expression (|8.14|l with the Schrodinger equation H8.13|l and using the 1 /Nc expansion for the 
energy 

Eb^N,Eo + J\Eb + 0{N-^), (8.16) 

we find in the leading order 

\ ^ V,nnpqTm.n{l)Tpq{l) = Eo , (8.17) 



2 

mnpq 



where 



In the next-to-leading order we obtain 



n- I 5Wbar(7) ,o 
Tmnil) = Z^lmj—^ ■ (8.18) 



E 

mnpq 



V ^ ( dluABii) , 1,, >^ dTpgi^) 

mnpq mn\i } / ^ "^pk r>. ^~ r) *^nmpq / ^ '^mj 



= ^Eb ■ (8.19) 



D. Next-to- leading order 

In the NLO we have equation H8.19|l . The structure of this equation allows for the solutions of the form 

~AB{l)=A^'\l)Yi\Ul)T- . (8.20) 
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AEb ^ AEq + UaAEa 

a 

where ^^"^^(7) obeys the inhomogeneous equation (|8.19|1 



E 

7nnpq 



mnpq mn\! ) / ^ ipk o ' n * ninpq / ^ fmj o 



and ^^(7) are sohitions of the homogeneous equation 

^ ^ ^mnpq'^mn (t) 
mnpq 



E 

L k 



Ipk- 



djqk 



^AE^. 



.21) 



.22) 



.23) 



E. Overlap matrix elements 



Now let us consider two baryon states \Bi) and \B2). Let us calculate the overlap matrix element for these two 

states 



^2■■■^N^ 



Ikjk 



/ d-id-i* exp -iVc E 1^3lh 



^ J d^dj* exp <-NAE -f^,% + [M^bar(7)]* + W^bar(7) ^B, (7) ^B. (7) 

(27Vc) I I J J 



/ d-fd-y* exp -Nc J2 Itjlt 



(8.24) 



At large Nc the integral over 7 can be computed using the saddle method. The saddle point equation is 

aM^bar(7) 



27* 



We use notation 7 = 7'"' for the solution of this equation. Combining Eqs. (|8.25|) and H8.15|l . we find 

y-^(o) (0)* ^ ^ 



(8.25) 



(8.26) 



F. From the equation for VFbar(7) to the Hartree equation 



Eq. I|8.17|l for W^bar(7) is vahd for any 7. If we take this equation at the saddle point 7*^"^ 
the equation 



, then we arrive at 



mnpq \ j 



(0) (0)* 
l^TpkTqk 



.27) 



This is nothing else but the expression for the Hartree energy (|2.14l) corresponding to the case M — 2, Lmn — 0: 



1 

'0 = 9 E E ^"1' 



r,s— 1 nin2n3n4 



(8.28) 
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Indeed, taking 



we find 



li^l = 4| E (<^™<)* (^12 = -£21 = 1), (8.29) 

r,s— 1 



E7S7S^*-^EC<- (8.30) 



Inserting this expression into Eq. H8.27|l . we reproduce the Hartree expression for the energy l|8.28|l . 



G. Saddle point equation for the baryon wave function 

We have explained above how to solve the Schrodinger equation for the baryon state \B) 1)8.1(1 described by the 
generating function <&b(7) (|8.5(I . Now let us turn to the vacuum. In our model the physical vacuum |0) state contains 
Nc quarks put into the bare vacuum \ This means that the physical vacuum of this model is described by the 
same equations as the baryon in simpler models with the trivial vacuum which were studied earlier. Therefore we can 
apply the old baryon representation ((3.3(1 to our new vacuum: 



|0)= Yl ^vac-'"=n<c|f^) 
ii...ijVc c— 1 

'l'vac(fc) — ^ ^ "^vac " ki-^ . . . ki^^ . 



ll...ZjV^ 



Here we use notation k for the "source" argument of the function $vac instead of g as it was in Eq. ((3.3(1 . 
Now we want to compute the matrix element 



.31) 
..32) 



n ^djajc \B) = 2^" ^ V'vac''"'=ffii (-Sbar)7i; 



No 



_d d_ 

dk* 



N, 



[$vac(fc)]*$s(7) 

fc=7=0 

[*vac(A:)]*$B(7) 



N. 



dQdQ*cxp I -^E'5y<3y 



El ^* d 1^9 



fc=7=0 



dQdQ* cxp 1^--^ E QijQ 



2 A^-^'^-^y 



fc=7=0 



Here 



N. 



Na 



JdQdQ*exp ^fj2QrjQh [*vac(g.g*)]*$B(Q) 



2 /^--cij-^ij 

1-3 , 



J dQdQ* cxp -^Y.Q^lQ 



2 ^ 



$vac(fc) = N^-A^^4k) exp [A^cWvac(fc)] 
^b(7) = ^r^B(7) exp [iV,I^bar(7)] • 



.33) 



(8.34) 
(8.35) 
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Note that the universal functional Wbar(7) describes the large Nc behavior of $5(7) for all low- lying baryons B. In 
the same way the functional Wvac(^) describes not only the exponential behavior of the vacuum $vac(fc) but also the 
"meson" functionals <I>incs(7) corresponding to the 0{N^) excitations above the vacuum. 
Thus at large Nc we have 



aNa 



:iV, 



J dQdQ* [Avac (g.g*)] AsiQ) expiVe <! [W^.dQg*)]* + W^UQ) - | E Q^iQl 
I »J 

JdQdQ*e^vl~fj:Q,,Q*A 



(8.36) 



At large Nc the integral over Q, Q* can be computed using the saddle point method. The deformation of the integration 
contour in the saddle point method can lead to the violation of the complex conjugation relation between Q^j and 
Q*y Therefore we introduce an independent notation for Q*^ 

8.37) 



The saddle point equations are 



dui 



[Wvac("*ff*)]* + W^bar(Q) " J I] Q 



mn '-^mn 



= 0, 



[T^vac(u*g*)]*+M^bar(Q)- 



Now we obtain from Eqs. 

These equations determine 

Combining Eqs. H8.42|l and (|8.15|) . we find 



2 ^ ^ Qmn^mn \ 
mn ) 



7 ^ * 



= 0. 



r5Wvac(fc)l 




rawvac(fc)i 


dki 







aVt^bar(Q) ^ 
^Q^J 

ij 



8.38) 

8.39) 

8.40) 
8.41) 

8.42) 

8.43) 
8.44) 



H. Functionals W{g) and Asig) from the saddle point method 
Applying the saddle-point method to Eq. H8.36|) . we find 

(0| fl {Y.9,a.j)j \B) = N:;-J{g) [A_Kg*)]* Ab{Q) 

X e^^Nc{[W.Uu*g*)r + W^Mi)}Q^Q(,)^^^^(,) ■ (8.45) 

Here J{g) stands from the contribution of the Jacobian corresponding to the fluctuations around the saddle point and 
Nc'^ accumulates all powers of Nc coming from various sources. We can rewrite the above representation in the form 

<0| it \yi9ja,c\ \B) = N:-AB{g) exp [NcW{g)] , (8.46) 



37 



where 

W{9) = + M/bar(Q.,)}Q=Q(,),„=„(,) , (8.47) 

Ab{9) = J{g) [Avac ^ ^ ^ . (8.48) 

Using the factorization of AsiQij) (|8.20(l . we arrive at the factorized form of Asig) 

Asig) ^ A^^\g)Y[Ak{g) (8.49) 

fc 

with 

Ak{g) = ^k{Q{g)) , (8.50) 

^W(5) = J(ff)i(°) (Q(5)) {Aac ([^^(5)5]*)}* • (8.51) 

Our result (|8.49|l shows the mechanism of the general factorization H1.18|l of functionals Asig) in models with the 
nontrivial vacuum. 

I. Representations for W(g) in terms of trajectories 

In our current model, the physical vacuum state contains Nc quarks occupying the bare vacuum. This means our 
physical vacuum can be described by equations derived in Sec. IVII for baryons. Thus the functional Wvac is given by 
Eqs. H6.39|l and l|6.42|l . In the special case Lmn = with which we are dealing now [See Eq. (|8.11(l ]. the result for 
Wvac is especially simple and is given by Eq. (|6.57|) : 

Wvac (k) = In /vac (fc) (if = 0) , (8.52) 
where /vac (5) is determined by the trajectories obeying asymptotic conditions H6.58|l . (|5.59|l and H6.60|l 

qTit) S,,oI..c{k) exp (slj) , (8.53) 

PTit) exp i-elj) , (8.54) 

9r(0) = fc„. (8.55) 

These trajectories obey condition (|6.25|l : 

Epr(i)'?r(i) = i- (8.56) 

n 

According to Eqs. H6.30|l . H6.41|) and 16.4211 we have 

pr(0) = . (8.57) 



One can generalize the work done in Sec. IVII for the case of states containing Nc quarks above the bare vacuum 
and derive a similar representation for Wbar(5): 

Wbar(g) =ln/bar(g) (if L„„, = 0) . (8.58) 

Here Qij is the boundary t = value of some trajectory Q^-"^{t) in the space of antisymmetric tensors 

Q^" (0) - Q^J ■ (8.59) 
and the parameter /bai (Q) is determined by the t — s- — oo asymptotic behavior of these trajectories 

gb-(i) '-^^ £,,/bar(Q) exp [(ei„ + eLr) t] , (8-60) 

P^'it) '-=°^ exp [- (eL.. + s^,) t] . (8.61) 
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Here Sij is the antisymmetric tensor 



1 ifi = l,j = 2, 
otherwise . 



.62) 



We assume that the values i = 1,2 correspond to the two occupied single-particle levels of the baryon state in the 
Hartree picture. The Hamilton equation describing the dynamics in the space of antisymmetric tensors Q^j^' , P}j^'^ 
will be discussed in Sec. |3 By analogy with Eq. (|8.57l) we have 



p,r(o) 



2 



r(Q) 



The analog of Eq. 



IS 



.63) 



.64) 



The constant on the RHS of Eq. (|8.63|) is fixed by relations (|8.15|) . 

The above representations (|8.52|) for Wvac{k) and H8.58|l for Whai{Q) are independent of each other. However, in 
the expression H8.47|l for the functional W{g) 



WsM = {Wvac [k{9*)]}* + W^bar [Q(g)] 



(8.65) 



the variables Q and k on the RHS become g (g*) dependent. We use the label Sch in Wsch(5) in order to stress 
that we are dealing with the Schrodinger approach. Later we will show how the functional W{g) can be calculated 
in the path integral approach [using notation Wpi(.g)]. Certainly the result of the calculation should not depend on 
the method but in order to avoid confusion at the intermediate steps we prefer to use different notations for the two 
different methods. 

The dependence of k and Q on g is described by equations H8.4U|I - H8.42|I : 



9W^ba.(Q) 





■9W^vac(fc)' 




■awvac(fc)" 


9j 


[ dk, \ 


- 91 


. 9kj _ 



Q^ 



(8.66) 
(8.67) 



Thus the problem of the calculation of the functional Wsch{g) reduces (in the case Lmn = 0) to search for trajectories 
obeying boundary conditions H8.53|l - H8.55|l . H8.59|l - H8.61|) . (|8.66|) . (|8.67|) . The asymptotic behavior of these trajectories 
at t ^ oo determines parameters Ivac{k) and /bar(Q)- According to Eqs. 18.52|l . (|8.58() and H8.65|l we have 



WsM - {In/vac [k{g*)]}* + In/bar [Q{g)] (if i™„ = 0) 



(8.68) 



J. Generalization for models with M > 1 



The previous analysis was devoted to "baryon" states containing 2Nc quarks above the bare vacuum. These states 
were described by functions ^sil) depending on antisymmetric tensors jij. The generalization to the case of states 
with MNc quarks is straightforward. In this case we must work with wave functions depending on antisymmetric 
tensors j[ij^i2...iM] rank k. The analog of Eq. (|8.8|) is 

^"1 d 

Tmn — / ^ O'mc^ric ~» 777 7TT Imi^ii-.-iM J (8.69) 



where the derivative is normalized by the condition 

d 



^nhh-jM = det||(5,„jj| . (8.70) 



39 



IX. PATH INTEGRAL APPROACH 
A. Advantage of the path integral approach 

Our previous analysis of the large- A^c models was based on the operator approach with the stationary Schrodinger 
equation (|2.2|l as a starting point. Now we want to study large- iVc models using the path integral method. In principle, 
both Schrodinger and path integral approaches must lead to the same results. In Sec. El we will demonstrate the 
equivalence of the results derived using these two methods. However, in spite of this equivalence we will see that 
the path integral method allows us to obtain in a straightforward way rather interesting results which are not so 
obvious in the Schrodinger approach. The reason is that the path integral approach is based on time-dependent 
trajectories. On the other hand, as we know from Sec. IVII in the Schrodinger approach the trajectories appear rather 
indirectly. In the Schrodinger approach one first derives the Hamilton-Jacobi equation for the "action" W{g) and 
only after that W{g) can be interpreted in terms of trajectories. In the case of models with the nontrivial vacuum 
the situation is even more involved, since the trajectories appearing in the Schrodinger approach belong to the phase 
space of antisymmetric tensors. On the contrary, the classical dynamics coming from the path integral approach can 
be formulated in simpler terms. 

Since the path integral approach allows for a straightforward derivation of the representation for the functional W{g) 
in terms of simple classical dynamics, we prefer to turn to the path integral formalism in this section. The equivalence 
of the results obtained in this section with the results based on the Schrodinger equation will be established in Sec. 

m 



B. Model 



We want to apply the path integral approach to models described by the Hamiltonian H2.8|l . For simplicity we will 
consider the case when L„j^„2 = in this Hamiltonian. For our aims it is convenient to rewrite the Hamiltonian in 
the form 

^ = ^ E («^r„a) Vc^p (a+Tpa) , (9.1) 



Q/3 



where 



(a+r„a)=^^a+„rr 



c— 1 myi 

The coefficients Vap and F™" are assumed to obey the conditions 

Vafj — V*p , Vap — V/3a ■ 
\^ a } ^ a 

The expression (19.11) for the Hamiltonian corresponds to the choice 

T/ — \ ^ -rnin2pn3n4 



(9.2) 



(9.3) 
(9.4) 



(9.5) 



in Eq. ((231) • 

In the path integral approach this system is described by the action 



dt 



(a+dto) 



2Nc 



af3 



(a+Fafl) Vap (a+F^a) 



(9.6) 



We use the same notation a"*", a for the Grassmann integration variables as for the operators in Eq. H9.1II . Strictly 
speaking, the transition from the operator formalism to the path integral has to be accompanied by a careful treatment 
of the problems of the operator ordering. However, in the leading order of the 1/Nc expansion these subtleties are 
not important. 
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One can "bosonize" the theory introducing the path integral over the auxihary boson variable Tr{t): 

J D7rcxp{- [a+K{TT)a]+NcShosiTT)} 



exp 5£;(a) 



Here 



/ DTrexp [Arc5'bos(7r)] 



Q/3 



Q/3 



[a+X(7r)a] = ^ dta+K{TT)a , 
K{7t) ^dt+J2 '^"^a . 



(9.7) 

(9.8) 
(9.9) 

(9.10) 
(9.11) 



C. Path integral approach at finite Nc 

The generating functional $ b (ff ) for the baryon wave function H2.6|l can be represented in the form 

*B(5) = (0|J(5,a,0)|i?), 

where we use the short notation 



J{g,a,ti) = Y[ 



c=l 



m 



We work with the Euclidean time. Therefore operators a(t) and a'^(t) are not Hermitean conjugate. 
In order to study ^sig) in the path integral approach we introduce the correlation function 

Zig, g', - (0|T {j{g, a, h)J{g', a+ M)] |0) . 

The path integral representation for this correlation function is 

/ DaDa+DnJig, a, h)J{g', a+ , h) exp {- [a+K{^)a] + N,S^,,{^)} 

^' /Dai?a+i?7rexp{-[a+if(7r)a]+7V,5bos(vr)} 

At large ti — ^2 ^ +oo only the contribution of the lightest baryon \B) survives in Eq. H9.15|l : 



Z{g,g',ti,t2) 



(0| J(.9, a, 0)|B)e-(*^-*^)(^--^™^)(i?| J(g', a+, 0)|0) 



(9.12) 

(9.13) 
(9.14) 

(9.15) 
(9.16) 

(9.17) 



D. Large Nc limit 

Now we want to consider the limit of large Nc- In this case we have with the exponential accuracy 

(0|J(5,a,0)|B)~e^=^(^). (9.18) 
The difference £b ~ Svac appearing in Eq. H9.17|l has the order 0{Nc): 

£b - fvac = {E^,, ~ £;vac) + 0{N°) . (9.19) 
If we combine the large-time limit with the large- A^c limit 

ti - ^2 -> oo, Nc^oo, (9.20) 
then we find from Eqs. ftTTTjl . (pO^ and I^J^ 

Z{g,g',h,h) - exp {N, [- {Ey,,, - E^^) {h - ^2) + W{g) + [W{g'*)]*] } . (9.21) 
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E. Calculation of the path integral 

First we can compute the Gaussian integral over quarks in Eq. (|9.16() 

/Di{[DetA»lexpSi..(ir)}"- 

At large Nc the path integrals in the numerator and in the denominator can be computed using the saddle point 
method. The numerator and the denominator have different saddle points which will be denoted tt"^^ and tt™'^, 
respectively. The result of the saddle-point integration is 

-^lnZ(g,5',ti,t2) = 5„onioc(7r=') + A5bos(7r^') , (9.23) 

where 

K{Tr) 



5no„ioc(^) = In I [3 • {ti\K-\TT)\t2) ■ g] Det-^ 



)i ' 



(9.24) 



A^bos (tt) = S^M - 5bos(7r™^) - / dt [Lbos(vr) - ibos(^™^)] • (9.25) 
The saddle point n'^^ is given by the equation 

= 0. (9.26) 



■7- [S'nonloc (tt) + AS'bos (t) 
OTT 



The solution tt'^' is t dependent. It also depends on g, g' ,ti,t2. 

n''^TT-\t\g,g',t,,t2). (9.27) 
On the contrary, the saddle point tt™'^ is a constant which can be found from the saddle point equation 

S 



. ln{[Detif(^)][exp5bos(^)]} 
ovr 



0. (9.28) 



7r— TT^ 



At large times i ^ ti or t <C i2 we have 



Combining Eqs. I9.21|l and (|9.23|) . we find 



lim TT^'it) ^tt"'"' . (9.29) 

i— »-±oo 



W{g) + [Wig'*)]* ^ lim [5nonloc(7r^') + A^bos (^^0 + (^1 - ^2) (^^bar - ^^vac)] ■ (9.30) 

F. Fermion formulation of the effective large- A^c theory 

In principle, Eq. H9.30|) combined with the saddle point equation H9.25|l contains everything what is needed for 
the calculation of the functional W{g). However, these equations involve the nonlocal action 5'nonioc(7r'^') defined 
in Eq. (|9.24l) . In order to solve the saddle point equation H9.26|l we would like to find a local formulation for our 
effective large- A^'c theory. To this aim we write the path integral representation for S'nonioc(7r) over auxiliary Grassmann 
variables b, 6+: 



exp^no„ioc(^) - [5 • {ti\K-\7r^^)\t2) ■ g'] Det- 

(9.31) 



A:(7r™'^) 

/ DbDb+ [gb{h)] [g'b+{t2)] exp S{b, b+, n) 



J DbDb+expS{b, 6+,7r™<=) 
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Note that the fermion fields 6, 6"*" appearing here have no color indices. The action S is 

5(5,6+,7r) = J dt[b+K{7r)b] . (9.32) 

The path integral lP0T|l can be also interpreted in terms of the operator approach to the theory of quantum fermions 
in the background classical field tt 

/i.W[,5(.0][,-5n^.)]exp5(M^^) ^ [ ] 3 

/D6i:)6+exp^(6,5+,^™=) \ i v lyj v i^j i v / 

For brevity we omit the in and ouf labels of the two "vacuum" states in this matrix element. One should also 
keep in mind that the Euclidean evolution violates the Hermitian conjugation of operator b{t) and b'^{t). Moreover 
even the hermiticity of the Hamiltonian can be violated. Indeed, the background field tt^^^ comes from the saddle 
point equation (|9.26() and the application of the saddle point method can be accompanied by the deformation of the 
integration contour. 

Combining Eqs. (|OT|l and ^^1^ . we find 

^„onioc(7r) = ln(0| [gb{h)] [g'b+it^)] |0), . (9.34) 

Inserting Eq. into Eq. we obtain 

W{g) + [W{g'*)]* = lim {ln(0| [gb{h)] [g'b+{t2)] |0),ci + A^bos (^^') - {h ~ ta) (^^bar ~ ^^vac)} • (9.35) 

Using Eqs. H9.8|) . H9.9|l . and (|9.34() . we find from the saddle point equation (|9.26() 
6 



(0| [gb{h)] [g'b+it2)] |0). \ = - E (^'% Wb^it^)] |0).c, . (9.36) 



The variational derivative on the LHS has a simple interpretation in terms of matrix elements in the operator formu- 
lation of the same theory: 

^^(0| [gb{h)] [g'b+it^)] |0). = -(0\T{[gbih)] [6+(t)r„6(t)] [g'b+it^)]} \0)^ . (9.37) 
Combining equations (|9.36f) and H9.37|l . we find 

E (^'% l9bih)] [g'b+ih)] |0).c, 

= (0|r { [gbih)] [b+{t)T^b{t)] [g'b+{t,)] } |0).c, . (9.38) 
This is a new form of the saddle point equation which determines tt'^'. 



G. Separation of the ti and t2 contributions 

Our aim is to compute the functional W{g). However, representation H9.35|l is written for the sum W{g) + [W{g'*)]* . 
Now we want to extract information about W{g) from this equation. At finite ti, t2 the RHS of Eq. (|9.35|) is a 
nontrivial functional of g and g' . Only in the limit ii — t2 — * +oo the dependence on g and g' reduces to the sum of 
two independent terms W{g) + [W{g'*)]* . This additive dependence is a consequence of the properties of the classical 
solution TT'^^{t) at large ti and t2- 

Eq. (|9.35|) was derived in the limit ii — t2 — +oo. But if we look at the history of the derivation of this equation, 
then we find that the limit ti — 12 ^ +oo was used only on the RHS of this equation. In particular, the saddle point 
equation l|9.38|l does not know anything about the limit ti — t2 — > +oo. 

According to Eq. H9.29|l at large t ±oo (when t ^ ti or t ^ t2) the field 7r'^'(i) approaches the constant 
configuration tt™'^. At large ti — t2 +oo we can also consider the intermediate region ±2 ^ t <^ ti. In this region 
the field 7r'^'(t) is also asymptotically static but with another value tt'^'*''. Thus we have three regions where 7r'^'(t) 
becomes asymptotically static: 

r 7r™^ ift«t2, 

7r'='(i) = <^ tt'^'^'- if ^2 < t < il, (9.39) 
[ 7r™'= ifti < 



43 



The constant configuration tt'"^'' corresponds to the "baryon" state \B) of the effective fermion theory (remember that 
now we deal with the effective fermion theory where fermions have no color indices). Therefore the matrix elements 
appearing in Eq. I)9.38|l become 



||r {[55(^1)] [g'6+(t2)]}|0).c 



[gbih)] \B).{B\ [g'b+ih]] |0).o, , 



(9.40) 



{0\T{[gb{h)] [b+{t)TMt)] [g'b+{h)]}\0)^.. 

{0\T {[gb{t,)] [b+{t)r^b{t)]}\B)^..{B\ [g'b+{h)] |0)^. . 
Inserting these factorized expressions into the saddle point equation H9.38(l , we find in the region ti ~ t 3> t2 



'\T{[gbih)] [b+it)T^bit)]}\B)^. 



(9.41) 



(9.42) 



Let us take the limit t2 —00 keeping ti and t fixed. In this limit Eq. (|9.42(l "forgets" about the behavior of 7r(t) 
at t ~ ^2 and at i <C ^2- Hence we can replace tt'^^ by its ti component 



As a result, Eq. H9.42|l becomes 



E ^9b{h)] - mmh)] [b+{t)r^b{t)] } 

/3 



Hence 



/3 



'\T{[gb{h)] [b+{t)rpbm\B)^a^ 



(0| [g6(ti)] 

Working with the background field tt'^^\ we can consider the limits t ±00 at fixed ti 

{0\T{[gbih)] [&+(t)r„6(0]} (0| [gb{h)] [6+(Or„6(0] 1^).- 



(9.43) 



(9.44) 



(9.45) 



(9.46) 



'\T{[gbih)] [6+(i)r„6(t)]}|i?),( 



Inserting these asymptotic expressions into Eq. 



[b+{t)TMt)] \0)r 
we obtain 



/3 



[gbih)] \B) 



(9.47) 



(9.48) 



;r^- = ^14^(51 [6+ (0)r^&(0)] \B)^. 

/3 



(9.49) 



In Sec. IIX HI we will show that these are nothing else but the static Hartrce equations for the vacuum and for the 
baryon in the large- iVc theory. 

Note that the time-ordered product has a discontinuity aX t = ti: 



lim — lim 

t— >ti-0 t^ti+O 



'\T{{9bih)] [b+it)r^bit)]}\B),a^ = (0| [gT^b{h)] . 
Combining this with Eq. H9.45|l . we conclude that the field n'^^^t) has a discontinuity at t ^ ti: 



lim — lim 



V-,. (0| [gr^6(ti)] 



(9.50) 



(9.51) 
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Using the property H9.40|l . we can rewrite Eq. (|9.35|) in the form 

W{9) + [W{g'*)\* = hi(0| [gb{h)] |B),<i, + \n{B\ [5'6+(t2)] |0),(2) + A^bos (^^') - {h - i2) (^^bar - ^^vac) ■ (9.52) 
Here -k^"^^ is the analog of tt'^^^ for the vicinity of the t2 time epoch: 



r 7r™= ift<t2, 
7r(2)(t) = <^ 7r^i(i) if t^t2, 
[ TT^^" ifi><2. 

Using the definition (|9.25() of AS'bos and the properties (|9.43() . H9.53() of tt^^^ and 7r^^\ we obtain 
A^bos (tt^') = 

J —OQ 

+ / dt [ibo.(7r^^) - ibo.(^™^)] + [ibos(^'^") - ibos(7r™^)] (ii - i2) 



[Lbos(7r^') - ibos(^™^)] + /*' dt [ibos(^^') - ibos(7r'^"'')] 



dt 



' — OO 

+ /" dt 

-/ — OO 



Lbos(7r(2))-Lbos(7r™^) 



+ 00 



dt 



dt 



Lbos(^('^)-ibos(^'^"^) 

+ [Lbo.(^'^"^) - ibos(7r™=)] (tl - t2) . 

Inserting this into Eq. 19.5211 and separating the ti and t2 contributions, we find 



Lbos(^('))-ibos(^'^"^) 

Lbos(7r(i))-Lbos(vr™^) 



T^(.9) = hi(0| [gb{h)] \B),^^, 



dt 



ibos(7r(i))-ibos(7r'^") 



dt 



ibos(^('))-ibos(^™^)' 



+ tl {(ii;bar - i?vac) + [ibos(^'^'^') " ibos(7r™^)] } • 



(9.53) 



(9.54) 



(9.55) 



Thus we have reduced the problem of the calculation of the functional W{g) to the analysis of the fermion system in 
the "self-consistent" field 7r(^^(t). Indeed, the field 7r'^^^(t) is determined by Eq. (|9.45|) combined with the discontinuity 
condition (|9.51() at ti and with the boundary conditions H9.43|l 



lim 7r(i)(t) = vr*^" , lim Tr^^\t) = tt" 

i— > — 00 t — )-+00 



(9.56) 



Inserting this field n'-^^t) into Eq. H9.55|l . we can find W{g). Note that the g dependence appears via the discontinuity 
(l^3T|) of 7r(i)(t) at tl. 



H. Hartree equations from the path integral formalism 

The saddle point field 7r'^^^(t) is determined by equation (|9.45|) . This equation contains the matrix elements 

(0| [56(ti)] , (9.57) 



{0\T {[gb{t,)] [6+(t)r^6(t)]}|B),a, 



(9.58) 



of the effective theory with the b fermions that have no color. 

The matrix elements (|9.57ll and (|9.58|l are written for the theory in the background t-dependent field tt'-^K The 
vacuum state (0| appears in these matrix elements as an out state corresponding to t — s- +00. At t — s- +00 we have 
7r'^^)(t) — !■ TT™'^ according to Eq. (|9.5()l) . Therefore the state (0| is the ground state of the theory with the Hamiltonian 
associated with the quadratic form l|9.11|l taken in the static background field tt^-^^: 



(9.59) 
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This means that the physical vacuum state (0| could be thought of as made of AI quarks above the bare vacuum {Q\ 



M 



M 



where the "single-particle" wave functions (f)^^^ are solutions of the equation 

a 

Using expression (|9.fi()|l , we can compute the matrix element on the RHS of Eq 

M 

{o\[b+{o)T^b{o)] |o).™. = ^(0:^j+r 

we find 



(9.60) 



(9.61) 



vac ■ 



(9.62) 



Inserting this result into Eq. 



M 



(9.63) 



Now we see that Eqs. ()9.61f) and (|9.63|l are nothing else but the Hartree equations H2.10|l . H2.11|l for the Hamiltonian 
(|9.1() . In the framework of these Hartree equations we have 



1 

-Bvac = ^ X! 



2 ^/ ^ ^vac ■ 

s-1 



(9.64) 



Instead of the second quantized representation for the physical vacuum |0) we can use the language the of A/-particle 
wave function written in terms of the Slater determinant 



|0) 



1 



/Ml 



1 det ||(0:,J I 

l<s,i<M * 



(9.65) 



Similar equations can be written for the baryon state \B) which is associated in the matrix elements (|9.57|l . ()9.58f) 
with the limit t — s- +oo: 



M+l 



'M+1 



\B) = n ('/'bar&+) i^^> = n 



(9.66) 



(9.67) 



M+l 



/3 s=l 



bar ' 



(9.68) 



M+l 



1 



(9.69) 



The state \B) corresponds to the Slater determinant 

1 



\B) —^= det ||(0bar)„ 

V/(M+1)! ^<s.^<M+l ll^^bar^n. 



(9.70) 
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We work with the pure 4-ferniionic Hamiltonian (|9.1|l which has no quadratic piece. Therefore we have additional 
simphfications. Combining Eqs. (jHS, f^^ . ^^U^, we find 



M 



^vac 



Q/3 



and similarly for the baryon 



M+l 



a/3 



Now we can simplify expression H9.55|l 

Wig)=ln(0\[gb{h)]\B),a^ +h 



M 



'M+l 

/ J ^bar / ^ ^vac 



dt 



(9.71) 



(9.72) 



(9.73) 



I. Time dependence 

In the previous section we have shown how the static Hartree equations for the baryon and for the vacuum appear in 
the context of the calculation of the generating functional W{g) for the baryon wave function. These static equations 
correspond to the asymptotic hmit H9.56|l of the field iT^^'^(t) at t — > ±00. 

Now we turn to the calculation of the time-dependent matrix elements (I9.57f) . (|9.58|l . These matrix elements 
are written in the Heisenberg representation with fixed states (0| and \B) but time-dependent operators h{t) and 
h'^{t). However, we can rewrite these matrix elements in terms of the Schrodinger representation. The Schrodinger 
representation is based on the time-dependent wave functions corresponding to the static Heisenberg wave functions 
(0| (liO)?! and \B) (ESU): 



M 



s=l 



(9.74) 



M+l 



(9.75) 



We use notation 4'vacit) instead of [0vac (*)]''' because the Euclidean evolution is not unitary. 

The time dependence of the wave functions 0vac(O '^bar(^) controlled by the quadratic form K^tt^^^) (I9.11|) 
which defines the b, theory 1)9.31(1 . 1)9.32(1 . Obviously the role of the Hamiltonian in the quadratic form K{n^^'') is 
played by 



The corresponding evolution operator is 



U{ti,T2) '^'=^^ Tcxp 



(9.76) 



(9.77) 



This operator U{ti,T2) controls the evolution of wave functions (p^^^{t) and (j)fj^j.{t): 



(9.78) 



lim I 

T— !- + CX3 



J+[/(T,i) exp(r<^J 



(9.79) 
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T^ — OQ 



(9.80) 



On the RHS we have the time independent wave functions (f>^^^ and (p'^^^. which are determined by the Hartree equations 
(inSU, J^THa for 0^^^ and (HTHTll . J?0)HI for 0^^^. 

For the matrix element H9.57(l . the reduction to the Schrodinger representation can be done as foUows 



M 



' M+1 



(0|[55(ti)]|S),a, =(0,ti|(55)|i?,ii).(i) -(f^ljn [-^vacliO^ljl^^^ | 1^^) • (9.81) 

Let us define 



git)=gU{h,t), 



(9.82) 



using the evolution operator H9.77|l . Next, let us introduce a special notation for the set of M + 1 functions including 

^vac 

{t) and g{t) 



0;,,(t) ifl <r<M, 
g{t) ifr = M + l. 



Then we find from Eq. (|9.81|l 

(oi [gb{h)] = (f^i <! n y>:^c,gih)b 



'M+l 



r=l 



\n) ^ det 15^5(^1) 

l<r,s<A/+l 



where we have introduced the notation 

Inserting Eqs. (^7^ . and using Eq. (|^75|l . we find 

Drs(t)^ lim lim exp(rie;^^-r2£^^,)(0^^J+C/(ri,r2)(/)^^ 

ri — ^+00 r2 — > — oo 



We see that Drs{t) is f independent 



Drs{t)=0, 



since the t evolution (|9.79() . H9.82|l of ^vac,g(^) compensates the t evolution (|9.80() of (j)^^^{t). 
By analogy with Eq. H9.81|l we can write 

(0|r|[56(ti)] [6 (t)r/36(t)J||i?).(i, -| (0,t|(&+r^6)[g(t)&]|B,t)^a) ifi>ti. 



In the case t < ti this leads to 

(0|T{[56(ti)][6+(i)r^6(t)]}|B)^a, {n\ I J] [<^;ac,s(ii)& 



' A/+1 



Inserting Eqs. 



X J\/+l 
^ r.s— 1 



and 



(M+i ^ 
'^pb) n [<^bar(ii)&+] 



\n) 



into Eq. (|9.45|) . we obtain 

M+l 
/3 r,s=l 



(9.83) 



(9.84) 



(9.85) 



(9.86) 



(9.87) 



(9.88) 



(9.89) 



(9.90) 
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According to Eq. (|9.85() Drs is a local (in t) function of 4>vac{t), 0bar(O and g{t). Therefore Eq. H9.90|l provides us a 
local expression for na{t) via (/)vac(^)i ^bai(0 ^^'^ 



where 



A/+1 |. 

ni") (0vac,0bar,<?) = E ^"/^ ^ [' 



<^vac,or/30ba 



(9.91) 



(9.92) 



The 0vac,g is defined similarly to Eq. H9.83|l . 

It is easy to see that the above argument can be repeated for the case t > ti: 



Inserting Eq. 



^W(i) '>J^ nw (0_(i),0ba..(t),5(O 

into Eq. ^7^, we find 



(9.93) 



W(g) = In det A-s + h 

l<r,s<Af+l 



'M+1 



M 



dt 



ibos(^(')) - ibos(^'^")] + r dt [Lbos(7r(i)) - Uosin^n 



(9.94) 



According to Eq. H9.87|l Drs is independent of t. However, Drs depends on ti. This ti dependence is implicit. Indeed, 
if we shift ti in the original path integral representation. 



then this will lead to the corresponding shift of the saddle point tt^^^^ 

^(i)(i)^^(i)(^_Ar) 

and to the shift in the evolution operator 

C/(ri,T2)->[/(ri-AT,T2-AT). 
Therefore the matrix Drs ()9.8t)|l will change as follows: 

Drs^Drs + AT{e:,,-eU . 



(9.95) 
(9.96) 
(9.97) 
(9.98) 



We see that the AT shifts of ti (I9.95|) and Drs (|9.98|) compensate each other in Eq. (|9.94|) so that the functional 
W{g) is independent of ti as it should be. Using this ti independence, we can simplify our equations by choosing 



ii = 0. 



Then Eq. H9.94|l becomes 



W(g) = In det Dr 

l<r,s<Af+l 



dt 



Lbos(7r(i))-Lbos(7r'^"'') 



dt 



Lbos(7r(i))-Lbos(vr™=) 



(9.99) 



(9.100) 



In the rest of the calculation we will assume the choice (I9.99f) . 

Using the expression (|9.9() for L^os and the result H9.90|l for tt'^^^ at i < 0, we find 



a/3 



(9.101) 
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Here 



Thus 



M+1 



r,s— 1 



We see that the RHS is t independent. Taking the hmit t — oo and using Eq. H9.56|) . we obtain 



Af+l 



(9.102) 



(9.103) 



(9.104) 



Combining Eqs. ^JU^, p.l04|l . we find 



In a similar way one can show that 



As a result, Eq. H9.94|l can be simplified: 



dt 



dt 



Lbos(7r(i))-ibos(7r'=") 



ibos(^('))-ibos(^™^) 



WpA.ig) =li^^ det D^s- 

l<r,s< Af +1 



(9.105) 



(9.106) 



(9.107) 



We use the label p.i. in order to emphasize that this resuh is obtained in the path integral approach. 



J. Closed system of differential equations 

The U evolution (I^TTjl of functions ^^ac(0> '/'bar(0> 9{t) described by Eqs. 1^1^ . and lOl can be 

represented in terms of the differential equations 



d 7. 



.it) 



Q 



dt 

a 

Using expressions ()9.91() and (|9.93|l for 7ri^^ , we can rewrite these differential equations in the form 

^^<PvacV 



.(OEni*^ (^vac(O,0bar(O,5(i)) 



^'/'ba. W = -E^L*' (0vac(i),'^bar(O,5(i)) r„</.L,.(0, 

a 

^5(0 = g{t) E ni*^ (<^vac(t), </.bar(t), <?(t)) , 



(9.108) 
(9.109) 
(9.110) 

(9.111) 
(9.112) 
(9.113) 



where Ili^'' should be used for t > and lli for t < (we work with ti =0). This is a closed coupled system 
of differential equations for functions 4'va.cit)i4'hiir(t),g{t)- In addition we have the boundary conditions which follow 
from Eqs. ^^7^ . (EHOI), and j^^H^ 



'^vac(*) ^ exp (ie^aJ (0^aJ , 



(9.114) 
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'^bar(i) exp i^teU </>Lr , (9.115) 



5(0) = 5- (9.116) 

Note that the g dependence appears via the boundary condition (|9.116() . Solving the system of differential equations 
(|9.111|I - H9.113|I with the boundary conditions H9.114|) - (|9.11fci|) . inserting the solution into Eq. H9.117|) and using Eq. 
(|9.85|) . we can compute the functional Wp ^Xg) H9.117|) : 

Wp.dg) = In det te,,,,(0) • <^Lr(0)l . (9.117) 

l<r,s<i\/-fl L 

This construction solves the problem of the calculation of the functional Wp. 1.(17). However, one can wonder how this 
representation for VFp.i.(5) derived in the framework of the path integral formalism is related to our results obtained 
in the Schrodinger approach. This question is studied in the next section. 

X. EQUIVALENCE OF THE SCHRODINGER AND PATH INTEGRAL APPROACHES 

A. Problem of equivalence 

Our analysis of models with the nontrivial vacuum has resulted in two expressions for the generating functional 
W{g). One expression (|8.68|l was derived from the Schrodinger equation using the operator formalism. The other 
result H9. 11711 for W{g) was obtained in the path integral formalism. Certainly both approaches are equivalent and 
should lead to the same results. However, in order to see this equivalence explicitly some extra work is needed. In 
Sec. IX Bl we demonstrate the equivalence for the simple case of models with the trivial vacuum and after that turn 
to the much more interesting models with the nontrivial vacuum. 

B. Models with the trivial vacuum 

Let us start from the simplest case when the physical vacuum coincides with the bare one. Then the set of occupied 
vacuum levels is empty 

{0:,,(i)} = (10.1) 

and 

^vac = . (10.2) 

The baryon is described by the Hartree equations H9.67|l - H9.69|l corresponding to one occupied level: 

^ (^^"^r„) 0bar = Ebar-^bar , (10.3) 



V.P (0Lr/30bar) , (10.4) 
b 



1 

Ehai = 2'^bar . (10.5) 
Since the physical vacuum coincides with the bare one, we have 

(0,t| = (0|, (10.6) 



51 



(0|6 = 0. (10.7) 

Therefore the matrix element (I9.88|) vanishes at i > ti. According to Eq. (|9.99|) wc work with ti = 0. Thus equation 
results in 



(0|T{[56(0)] [b+{t)rpb{t)]}\B)^n^ - 0(-i)(O| (b+Tpb) \B,t)^a^ = 0{-t) [g{t)r0<l>^,,{t)] . (10.8) 

On the other hand, Eq. (|9.84() gives 

(0| [gb{0)] = 5(t)0bar(t) = const , (10.9) 

which is t independent. Now we insert Eqs. pu.8|l and (|10.9|) into Eq. I9.45|l : 

K W = 0[-t) > Vaf3 r , . T-vp. (10.10) 

^ [5(t)0bar(i)J 

According to Eqs. H9.109|) and (|9.11U|I we have the differential equations 

^<^bar(i) = -^^i')(i)ra0bar(i), (10.11) 



j^9{t)^9{t)Y.^^^\t)T^ (10.12) 

a 

with the boundary conditions l|9.115|) . H9.116|l 

0bar(i) exp (-tebai) 0bar , (10.13) 

g{0) = g. (10.14) 
According to Eqs. H9.56|l . (|10.3() . and (|10.12() the asymptotic behavior oi g{t) at t ^ —oo is 

g{t) *^=°° const exp (tebar) (10.15) 
and the path integral result (|9.117|l for W{g) becomes 

Wp.i.Cg) =ln[g(i)0bar(i)] . (10.16) 

Thus the path integral calculation of the functional Wp,i,{g) is reduced to solving differential equations pU.ll|) . (|10.12f) 
with TT^^^ given by Eq. pU.10|) and with boundary conditions I|1U.13|) . (|1U.14|I . (|10.15|) . 

According to Eq. pu.9|l the RHS of Eq. (|10.16|l is t independent. Taking the limit t ^ —oo and inserting the 
asymptotic behavior (|10.15() . we can determine the constant on the RHS of H10.15() : 

g{t) ' 0+ , exp (tebar) • (10.17) 



Let us show that this path integral result Wp. 1.(17) is equivalent to the Schrodinger approach discussed in Sec. I VII 
Note that our path integral treatment dealt with the case Lmn — 0. In the Schrodinger approach this case was 
discussed in Sec. IVI Fl 

In order to match the equations of Sec. IVI with the equations of the path integral approach one has to use the 
following "dictionary" connecting the two formalisms: 

qit)^git), (10.18) 

p{t) ^ e-'^-'-^^^Ub.At) , (10.19) 
£0 = Ebar , (10.20) 
Hg) = [9imb.r{t)] . (10.21) 

On the LHS we have the quantities appearing in Sec. IVI Fl within the Schrodinger approach whereas the RHS is 
represented by the objects which were used in the path integral method. 

Now it is easy to see that the boundary conditions Ht).58|) . Ht).59|l . and (|6.t)U|l correspond to Eqs. H1U.18|I . H1U.19|I . 
and H1U.14|I . respectively. The Hamilton equations Ht).5|l are mirrors of Eqs. H1U.11|I - H1U.12|) . 
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C. Variables Q^^'> , P^^' 



In the previous section we have demonstrated the equivalence of the Schrodinger and path integral approaches in 
the case of systems with the trivial vacuum. Now we want to consider the general case when the vacuum is nontrivial. 

Our main lesson from the analysis of the systems with the trivial vacuum is that the equivalence between the two 
approaches can be established on the basis of the Hamiltonian description of trajectories used for the calculation of 
the functional W{g). 

In the case of the nontrivial vacuum we already have the representation for W{g) in terms of trajectories obeying 
the differential equations (|9.1U8I) - (|9.11U|I . Now we want to rewrite these differential equations in the Hamiltonian 
form. 

To this aim we must introduce the canonical variables. Note that the original differential equations H9.1()8|l - H9.110() 
are written in terms of single-particle wave functions (f>^^^ (t) , (t) . As we will see below, the Hamiltonian formalism 
is based on M- and {M + l)-particle wave functions including the Slater determinants associated with the states (19 .741) 
and 

M+l 

PL'!^X,At) = Antisym J] i^LMm, = 777-^ det l^LM^. , (10.22) 



M 

Ql^ltmjt) = AntisymH feac(0' 



= TTT det 

m. Ml l<r,s<M 



(10.23) 



We assume that the antisymmetrization operator Antisym is normalized by the condition 

(Antisym)^ = Antisym . 



(10.24) 



Assuming the choice H9.99II ti = 0, we deal with the M-particle states at i > and with the {AI + l)-particle states 
at i < 0. Therefore the M+l particle "wave function" Pmi^.mM+i is relevant for the description of the region t < 
whereas the variable qI^^.thm "^i^ used at t > 0. 

Now we want to introduce the canonically conjugate variables for Qm[\-rnM ^^id for Pmi.tmM+i- 



M+l 



PLV..mJt) = iM+l) ^ g™,,^,(i) Antisym [] [cj^LMr 



mi...mM + l 



(10.25) 



= Antisym g^,,^,{t) J] feac(i) 

mi...mM + l I - '- ■ 



s=l 



(Af + 1)! 



■det 



'vac.g (0 



(10.26) 



On the RHS of the Eq. (|10.26(l we use notation ^^ac.gC*) JHSSll- 



D. Differential equations for Q'^', P'^' 



Now let us rewrite the differential equations H9.108I '1~ (|9.11U|I in terms of variables Q^^'^\ p(A^) for t > and in 

terms of Q(*^+i), p(A^+i) for t < 0: 



(10.27) 
(10.28) 



where 



L = 



M ifi>0, 
M+l if i < , 



(10.29) 
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and 



^ ] ^i^a)m-^n -^nm2...mL^^'^ ^ ^^<^'^m2n ^min...mL^^^ ' (10.30) 

n 

|^'9nm2...mi, (0 i^a)nmi QmiTi...mi (0 (ra)„,„2 + • • ■ + Qmi...mL-ini^) i^a)nmL ' (10.31) 



The field Tri^-* (t) is given by 



" ^^"^ [QW(t).pW(t)] 



where L is defined by Eq. pu.29|l and 

n(i) . p{L)] - V 0^^'> P^^'> 

^ ^ / , ^mi...mL mi.. 



Obviously 



and 



knm2...mL 



It V 



(M +1) ■ P^*'+''(0)j - [q(*^H0) • P(*^)(0) 



(10.32) 

(10.33) 
(10.34) 

(10.35) 
(10.36) 



mi...mM + l 



E. Hamiltonian interpretation of differential equations 

Let us introduce the Poisson bracket 



{^5mi...mz, 1 -fyii-.-nL } 7-1 ^Ct (5772-72^ 
L\ jk " 



Then 



{4^U,[Q(^^-r,-p(^)]} = -[r,.g 



7111 ...rriL 



Now Eqs. (|10.27() and H10.28|) take the form 



Qm}m2...mL ^ ] (0 |^Qmim2...m_L ' IQ ' ^ ' J' 

a 

'^^mim2...mL ^ ^ T^ait) ^Pmim2...mL ' ' ' ^1 } 
a 



Inserting 7ri^^ from Eq. (|10.32|l . we can rewrite these equations in the form 



^^^mim2...mL 



mim2...mi. 



Sm-im2...mL i ^ f 

[QW{t)-PW{t)\ 



[Q(^)(t)-PW(t)] 



(10.37) 

(10.38) 
(10.39) 

(10.40) 
(10.41) 



(10.42) 
(10.43) 
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where the Hamiltonian H is 



(10.44) 



Q/3 



According to Eq. (|10.35|) the denominator \Q''^\t) ■ P''^\t)\ is t independent along any solution so that this denom- 
inator can be treated as a constant. 



F. Boundary conditions 
Using Eq. 1)10. 2H|I . we can rewrite the boundary condition l|9. 114(1 in the form 



Similarly we find from Eqs. (|9.115|l and (|10.22|) 

p(M+l) 
^ mi...mM + i V / 



t — > — oo 



det ( 

<r,s<M+l 



bar/m. 



Comparing Eqs. (|10.23|) and pu.26|) . we find at t = 



Q^f.+'i.,+,(0)= Antisym [Qltf.l™., (0)5mM+i 

mi...mM+l '- 

Similarly, the comparison of Eqs. pU.22|) and pu.25|l results in 

niM + l 



(10.45) 



(10.46) 



(10.47) 



(10.48) 



According to Eqs. 



G. Expression for W{g) 
mrm . and we have 

det Drs = {M + 1)! [q(*^+i) (t) ■ (t) 

l<r,s<M+l L 



Combining this with Eq. H9.107|l . we find 

Wp.i.(ff) = ln{(M + 1)! [Q(*^+i)(i) • P(*^+I)(i)] } 



(10.49) 



(10.50) 



Later we will check the equivalence of this result with the expression H8.68|l obtained in the Schrodinger approach. 

Note that the RHS of Eq. (|10.50|l is t independent according to Eq. (|10.35|l . Using Eqs. H10.35|) and (|10.36|) . we 
can also write 



p.,(5) = ln{M! [Q(*^)(<).p(*^)(i)]} 



(10.51) 



H. Asymptotic representation 



The asymptotic behavior of P^Lmui'^) at t ^ +oo is given by the expression 



'^/(5)exp (^-tEeLc) ^^det^^l 



(10.52) 
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with some constant I{g). According to Eq. H10.35() the combinati on \Q''^ \t) ■ is t independent so that it 

can be computed via its Umit at t —* +00. Using Eqs. H1U.45|) and (|10.52|l . 



— hm 

t — *+oo 



Q(Af)(i).p(M)(^) 

Combining this with Eq. H10.36|l . we find 

g(A^+i)(i).p(Am)(^) 



we obtain 

Q(M)(i).p(M)(^)l ^j(^gy 



1(9) 



M + 1 

Taking the hniit t — s- — 00 in this equation, we can fix the constant in the asymptotic expression 



)(M+i) lig) 

imi...mM+l y I 



Inserting expression (|10.53|l into Eq. (|10.51|l . we obtain 

VFp.,(.9)=ln/(.9), 



(10.53) 



(10.54) 



(10.55) 



(10.56) 



I. Comparison of path integral and Schrodinger results 

We have computed the functional W{g) using the operator (the Schrodinger) approach and the path integral 
method. The two results are given by Eqs. (|8.()8|) and H1U.56|) . respectively. Our path integral result (|10.56|l is more 
general since it deals with models where the physical ground state contains N^M quarks with arbitrary M whereas 
the Schrodinger result H8.68|l corresponds to M — 1. 

Comparing equations (|8.t)8|) and H1U.56|I . we see that they give the same result for W{g) if 



/(ff) = {/vac (fc(g*))}*/bar {Q{g)) 



(10.57) 



The functional I{g) arises in the path integral approach whereas the functionals /vac (^(ff*)) a-nd /bar ((5((7))come from 
the Schrodinger approach. The aim of this section is to give an independent derivation of the identity (|10.57(l which 
directly establishes the equivalence of the path integral and Schrodinger methods. 

The functionals appearing in Eq. (|10.57|) can be represented in terms of trajectories. In the path integral approach 
the trajectory representation for I[g) arises via the saddle point method. The Schrodinger representation for W{g) 
was formulated in Sec. IVIII II in terms of trajectories (7™'^(i), pl^'^{t) for the vacuum sector and trajectories Q^j^'X^), 
P'°^^{t) for the baryon sector. Let us show that the connection between the trajectories of the Schrodinger approach 
and the trajectories P'^^\t), Q'^^^t) (with L = 1, 2) of the path integral approach is given by equations 



Qf\t) = i:,Apr{-t)] 



(10.58) 
(10.59) 



QlfW- 2^:acQ^"W, (10.60) 
p'^^it)^h^P^(t). (10.61) 

Remember that the trajectories qj^^{t), p^'^^{t) were introduced in Sec. IVIII II for t < whereas the trajectories P'^^^t), 
Q'^-'(t) are defined at i > 0. Note that the equations for the trajectories are invariant under the combination of the 
time reversal transformation t ^ ~t with the complex conjugation and with the exchange P^^^ <-> Q^^\ Therefore 
functions P'^^^(t), Q'-^-'(t) defined by Eqs. (|10.58|) and 110.59|l automatically obey the equations of motion. 
Next, the equation for trajectories are invariant under the rescaling transformations 

Q^^^^aLQ^^\ P^^^^hLP^^K (10.62) 

Therefore various t independent factors appearing in Eqs. H1U.58|) -( |10.61|I do not violate the equations of motion. 

Let us summarize. Eqs. H10.58|) - (|10.61|l express the trajectories Q^^\t), p(^'(t) of the path integral approach via 
the trajectories g™'^(t), p™'^(— i), Q^'^'' {t), (5'^"(t) of the Schrodinger approach. If the Schrodinger trajectories obey the 
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equations of motion associated with the Schrodinger approach then the trajectories Q^^\t), P'^^^t) defined by Eqs. 
(|10.58|I - H1U.61|I satisfy the classical equations motion I|1U.42I) . I1U.43II corresponding to the path integral approach. 

The next step is to check that the boundary conditions H1().45|I - H1().48(I of the path integral approach also follow 
from the equations of the Schrodinger approach. 

According to Eq. H8.55|l we have 



P^'^(O) -/barA:*. (10.63) 

Using Eqs. (|F37|) and ((TU3gjl . we find 

9Wvac(A:„) 1 



dkn /, 

Let us insert this equation into Eq. (|8.67() : 

1 



(0) . (10.64) 



gjQ'l^ (0) - g^Qf^ (0)J = . (10.65) 

According to Eqs. H8.59|l and (|10.61|) we have 

/4,Q., =2Qg)(0). (10.66) 

Thus 

5,Qf^(0)-g.Qf (0) = 2Qif(0). (10.67) 

This is nothing else but the boundary condition (|10.47() of the path integral approach. 

Now let us check the boundary condition H1U.48|I . We start from equation (|8.66|) and insert k from Eq. H1U.63|) 

E%^5. = y^/^^^^(0). (10.68) 



3 



Next we combine Eq. with Eq. (|10.61|l 



%^ = y^/^?(0). (10.69) 

Cf^ij -'bar 

Inserting Eq. (|1U.69|1 into Eq. (|1U.68|1 . we find 

2Y,pg\0)g,^Pl'\0). (10.70) 
j 

This shows that the boundary condition H10.48() of the path integral approach also holds. 
Combining Eqs. with Eqs. ((TIT3^ . ((TU3S|) . we find 

P.i'^ (t) * (5„o/bar/vac CXp {-el,,t) , (10.71) 

Qi'Ht)'^^°^ d,,oexp{el,t) . (10.72) 

The last equation agrees with Eq. (|10.45|l . 

Next, according to Eqs. (IHUl), H10.60|l . (|10.61|l we have 



Q[f (t) *^=°° ^e.jI'Lch.. exp [(eLr + eL) t] , (10.73) 



The last equation agrees with Eq. (|10.46|l . 

Thus we have checked that Eqs. (|10.58|) - (|10.61|l connecting the trajectories of the path integral approach with the 
trajectories of the Schrodinger approach are compatible both with the equations of motion and with the boundary 
conditions. 
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Now we can turn to the proof of the equahty of functionals Wsch{g) and Wp,i,{g) associated with the Schrodinger 
and path integral methods. As was explained above, this proof reduces to the derivation of the identity (|10.57|) . Using 
Eqs. CUSi, mmj, dUni we can compute I{g) XTrTh^. fTTT^ 



i{g) = [g(i)(i) • p(i)(i)J = 2 [g(2)(i) . p(2)(^)j ^ . 

This proves relation pu.57|) . Thus 

WsM = Wp.iig). 

The equivalence of the Schrodinger and path integral representations for W{g) is checked. 



(10.75) 



(10.76) 



J. Analyticity and time contour 

The functional $5(5) is holomorphic in g by construction (j2.6|l . This analyticity is inherited in the large- iVc limit 
(I2.7|l by the functional W{g), although this limit can be accompanied by the appearance of singularities in W{g). 
The representation for W{g) in terms of classical trajectories is compatible with the analyticity of W{g). Indeed, the 
Hamiltonian (|10.44(l used in this representations is a holomorphic functions of coordinates and momenta. However, 
the configuration of trajectories used for the construction of W{g) can lead to the singularities in W{g). Examples of 
these singularities can be found in the analytical expressions for W{g) in the rotator model analyzed in Sec. IVIII [see 
Eq. (|7.74|) ] and in the naive quark model studied in Ref. (16] . 

The analyticity of W{g) in g is closely related to the analyticity of the corresponding trajectories in t. In our 
path integral construction we used Euclidean time. However, the choice of the Euclidean time is not necessary. For 
example, in the Schrodinger approach the functional W{g) can be represented in a form which requires no time and 
no trajectories. The Schrodinger representation is formulated in terms of the action which is defined without explicit 
usage of the time variable. But if one wants to represent this action in terms of trajectories, then one has to use some 
version of the time variable. Naively the choice of the time contour (Euclidean, Minkowski or more complicated) is 
not important. However, even in simple systems we can meet singularities restricting the freedom of the choice of the 
time contour. Our general representation for the functional W{g) in terms of trajectories ignored this problem which 
requires an additional analysis in any particular model. 



XI. CONCLUSIONS 



A. Main results 



The main aim of this work was to check the consistency of the assumptions used for the construction of the 1 /Nc 
expansion for the baryon wave function in QCD. We have concentrated on those nonperturbative features which 
cannot be directly proved in QCD. On the other hand, these properties can be studied in detail in models preserving 
the basic structure of the 1/Nc expansion in QCD. Our check is quite successful: in the analysis of models we could 
explicitly demonstrate all basic features: 

- exponential large- A^c behavior of the generating functional for the baryon wave function, 

- universality of this behavior, 

- factorization of the preexponential terms. 

In the absence of the direct proof of these properties in QCD, this model analysis combined with the independent 
arguments of Refs. |l6lll7j | (based on the evolution equation, spin-flavor symmetry and soft-pion theorem) gives us a 
rather solid self-consistent picture of the baryon wave function in large- Ac QCD. 

The main criterion for our choice of models was the possibility to perform the direct calculation of the 1/Ac 
expansion. In the simplest models we could obtain analytical results. In more general models the problem has been 
reduced to the solution of a coupled system of differential equations corresponding to the large- Ac effective classical 
dynamics. 



B. Phenomenological applications 



Although the construction of realistic models was not the aim of this work, our analysis still has a direct relation 
to phenomenology. 
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1) The practical applications of the 1/iVc expansion which can be of interest for phenomenology are based on the 
large-A'^c contracted spin-flavor SU{2Nf) symmetry. In the problem of the baryon wave function we have a rather 
nonstandard manifestation of this symmetry ■ If we assume this realization of the spin- flavor symmetry, then we 
can derive interesting results like the connection between the distribution amplitudes of nucleon and A resonance. 
The aim of this paper was to justify this assumption by checking the structure of the 1/Nc expansion with a direct 
calculation. 

2) The complete analysis of the effects related to the spin-flavor symmetry requires a thorough study of the next- 
to-next-to leading order (NNLO) of the 1/Nc expansion. In the current analysis we did not go so far. Most of the 
work was concentrated on the leading order (LO). We could also check the factorization properties appearing in the 
next-to-leading order (NLO). However, the subtle NNLO contributions crucial for the understanding of the spin- flavor 
symmetry effects remained beyond the scope of this paper. Nevertheless the successful LO and NLO check of the 
picture suggested in Refs. [la, [13 gives us additional evidence for the validity of this picture. 

3) Although the results of this paper are presented in a simple quantum-mechanical form, they can be directly 
rewritten in field-theoretical terms relevant for the phenomenological models with four-fermionic interaction of the 
Nambu-Jona-Lasinio type. 



C. Other aspects 



The methods used in our analysis of large- A^c models are also interesting from the point of view of the relation 
between the large-A systems and classical dynamics. Our study of models has also revealed a connection between the 
equations for the generating functional W{g) and such traditional equations of the many-body physics as the Hartree 
equation (including its time-dependent version) and RPA equation. Although the methods developed here cannot be 
applied directly to large-A'c QCD, they can be used in some special QCD problems like the diagonalization of the 
matrix of anomalous dimensions of the leading twist baryon operators [iq . 
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